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1. Introduction 

In the Standard Model, SM, transition rates of semileptonic processes such as — > uHu, 
with {u^ ) being a generic down (up) quark, can be computed with high accuracy in terms 
of the Fermi coupling Gp and the elements Vji of the Cabibbo-Kobayashi Maskawa (CKM) 
matrix |j^] . Measurements of the transition rates provide therefore precise determinations 
of the fundamental SM couplings. 

A detailed analysis of semileptonic decays offers also the possibility to set stringent 
constraints on new physics scenarios. While within the SM all — > uHi^ transitions are 
ruled by the same CKM coupling Vji (satisfying the unitarity condition \ Vik\'^ = 1) and 
Gp is the same coupling appearing in the muon decay, this is not necessarily true beyond 
the SM. Setting bounds on the violations of CKM unitarity, violations of lepton universality, 
and deviations from the V — A structure, allows us to put significant constraints on various 
new-physics scenarios (or eventually find evidences of new physics). 

In the case of leptonic and semileptonic K decays these tests are particularly signifi- 
cant given the large amount of data recently collected by several experiments: BNL-E865, 
KLOE, KTeV, ISTRA-F, and NA48. These data allow to perform very stringent SM tests 
which are almost free from hadronic uncertainties (such as the fi/e universality ratio in 
K£2 decays). In addition, the high statistical precision and the detailed information on 
kinematical distributions have stimulated a substantial progress also on the theory side: 
most of the theory-dominated errors associated to hadronic form factors have recently been 
reduced below the 1% level. 

An illustration of the importance of semileptonic K decays in testing the SM is provided 
by the unitarity relation 

\Vud\^ + iKsP + iKfel' = 1 + CNP . (1.1) 
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Here the Vji are the CKM elements determined from the various processes, having 



fixed Gp from the muon hfe time: = 1.166371(6) x 10~ GeV~ ||2|. eNP parametrizes 
possible deviations from the SM induced by dimension-six operators, contributing either 
to the muon decay or to the — > transitions. By dimensional arguments we expect 
eNP ~ M^/A^p, where Anp is the effective scale of new physics. The present accuracy on 



\Vus\, which is the dominant source of error in (1.1), allows to set bounds on enp around 
0.1% or equivalently to set bounds on the new physics scale well above 1 TeV. 

In this note we report on progress in the verification of the relation (^]^) as well as 
on many other tests of the SM which can be performed with leptonic and semileptonic 
K decays. The note is organized as follows. The phenomenological framework needed to 
describe K£-^ and Kf^2 decays within and beyond the SM is briefiy reviewed in Section |2|. 
Section^ is dedicated to the combination of the experimental data. The results and the 
interpretation are presented in Section ^. 

2. Theoretical framework 

2.1 and rates within the SM 

Within the SM the photon-inclusive K^^ and K12 decay rates are conveniently decomposed 

as m 



miSi,)) = ^Ucx5ew|Ks|V+(0)2 4(A+,o) {l + S§ui2)+Se^Y > (2-1) 



Vu. 



ud 



Ik^k /I - mym 



^ > X (1 + 5em) , (2.2) 



where Ck = 1 (1/2) for the neutral (charged) kaon decays, /|j-(A+,o) is the phase space 
integral that depends on the (experimentally accessible) slopes of the form factors (generi- 
cally denoted by A+^o); and 5ew = 1.0232(3) is the universal short-distance electromagnetic 
correction computed in Ref. [Q] . The channel-dependent long-distance electromagnetic cor- 
rection factors are denoted by (5em and 6^^. In the case 6em = —0.0070(35) [^, |6|, while 
the four 6^ are given in Table ||, together with the isospin-breaking corrections due to 
rriu / ma, denoted by Sfjj^^)- 

The overall normalization of the K^^, rates depends upon /+(0), the K ^ tt vector 
form factor at zero momentum transfer [t = (px — Pn)'^ = 0]. By convention, /+(0) is 
defined for the — > vr" matrix element, in the limit = rrid and aem (keeping 
kaon and pion masses to their physical value). Similarly, fx/fn is the ratio of the kaon 
and pion decay constants defined in the = and Oem — *■ limit. The values of these 
hadronic parameters, which represent the dominant source of theoretical uncertainty, will 



be discussed in Sect. 4.2 



The errors for the Ki^ electromagnetic corrections, given in Table have been obtained 
within ChPT, estimating higher-order corrections by naive dimensional analysis [0, |^. 
Higher-order chiral corrections have a minor impact in the breaking of lepton universality. 
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'^5C/(2)(^o) 









+0.57(15) 




2.36(22) 


+0.08(15) 


K% 




2.36(22) 


+0.80(15) 
+0.05(15) 



Table 1: Summary of the isospin-breaking corrections factors |0, ||]. The electromagnetic correc- 
tions factors correspond to the fuUy-inclusive -fi'£3(^) rate. 



The errors are correlated as given below: 



/ 1.0 0.1 0.8 -0.1 \ 
1.0 -0.1 0.8 
1.0 0.1 
1.0 



(2.3) 



2.2 Parametrization of Ki^ form factors 

The hadronic IT — > vr matrix element of the vector current is described by two form factors 
(FFs), U{t) and/o(t), defined by 

(vr^ (k) |s7^n|i^°(p)) 
f-it) 



{v + krf+{t) + {p-kYf-{t) 



m 



(2.4) 



t 



where t = [p — k)^. By construction, /o(0) = /+(0). 

In order to compute the phase space integrals appearing in Eq. ( |2.1| ) we need ex- 
perimental or theoretical inputs about the i-dependence of f^^o{t). In principle, Chiral 
Perturbation Theory (ChPT) and Lattice QCD are useful tools to set theoretical con- 
straints. However, in practice the i-dependence of the FFs at present is better determined 
by measurements and by combining measurements and dispersion relations. 

In the physical region, (m| <t< {itik — ?n.7r)^), a very good approximation for the 
FFs is given by a Taylor expansion up to terms 

/+,o(t) 



/+(0) 



1+V ^ +1 A" r * 



-,0 2 

' mi 



+ 



(2.5) 



Note that t = {pK — P-k) — ""^i: + ~ 2mKET^, therefore the FFs depend only on E^^. 
The FF parameters can thus be obtained from a fit to the pion spectrum which is of the 
form (/(-Ett) X /(i?^)^. Unfortunately t is maximum for E.,^ = 0, where g{ET^) vanishes. 
Still, experimental information about the vector form factor /+ measured both from 
and K^^ data are quite accurate and so far superior to theoretical predictions. A 
pole parametrization, f+{t) = My /{My — t), with My ~ 892 MeV corresponding to the 
-Rr*(892) resonance and which predicts A" = 2(AY)^, is in good agreement with present data 
(see later). Improvements of this parametrization have been proposed in Refs. |^, |l^, |TT|l . 



For instance, in Ref. |11|, a dispersive parametrization for /+, which has good analytical 
and unitarity properties and a correct threshold behavior, has been built. 
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The situation for the scalar form factor /o(t) is more complex. For kinematical reasons 
fo{t) is only accessible from K^^ data and one has to deal with the correlations between 
the two form factors. Moreover, for foit), the curvature Aq cannot be determined from the 
data and different assumptions for the parametrization of /o such as linear, quadratic or 
polar lead to different results for the slope Ag which cannot be discriminated from the data 
alone. In turn, these ambiguities induce a systematic uncertainty for Vus, even though data 
for partial rates by itself are very accurate. For this reason, the parametrization used has 
to rely on theoretical arguments being as model-independent as possible and allowing to 
measure at least the slope and the curvature of the form factor. 

2.2.1 Dispersive constraints 

The vector and scalar form factors /+,o(*) i^i Eq. (|2.4| ) are analytic functions in the complex 
t-plane, except for a cut along the positive real axis, starting at the first physical threshold 
*th = {nT-K + "ivr)^) where they develop discontinuities. They are real for t < t^h- 

Cauchy's theorem implies that f+fl{t) can be written as a dispersive integral along the 
physical cut 

oo 

f,Q(t) = -[ ds' \ + subtractions , (2.6) 

TT J [s' — t — iO) 

where all possible on-shell intermediate states contribute to its imaginary part ImFfc(s'). A 
number of subtractions is needed to make the integral convergent. Particularly appealing 
is an improved dispersion relation recently proposed in Ref. |12] where two subtractions are 



performed at t = (where by definition, /o(0) = 1) and at the so-called Callan-Treiman 
point tcT = {f^K ~ "^V) leading to 



/o(t) = exp 
with G{t) 



t 

tcT 

tcT{tcT-t) rds' 



In (/o(tcT)) - Git)) 



(2.7) 



^ {s' -tcT){s' -t-ie)' 



assuming that fo{t) has no zero. Here (p{x), the phase of fo{t), can be identified in the 
elastic region with the S-wave, 7 = 1/2 Kir scattering phase, Sk-k{s), according to Watson 
theorem. 

A subtraction at tcT has been performed because the Callan-Treiman theorem implies 

where Act ~ (^{nT'u,d/^'^F'n) is a small quantity. ChPT estimates at NLO in the isospin 
limit obtain 

Act = (-3.5 ± 8) X 10"^ , (2.9) 

where the error is a conservative estimate of the high-order corrections to the expansion 
in light quark masses [^]. A complete two-loop evaluation of Act, consistent with this 
estimate, has been recently presented in Ref. [p 
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Hence, with only one parameter, foitcr)-, one can determine the shape of /o by fitting 



the i^^3 decay distribution with the dispersive representation of fo{t)^ Eq. (2^). Then, we 
can deduce from Eq. ( |2.7D the three first coefficients of the Taylor expansion, Eq. ( |2.5| ), see 
Ref. |T|: 



A 



K-n 



In /o(tcT) -G(0) 



In /o(tCT) - 0.0398(40)) 



A^' = (A^)^ - 2 ml/tcT G'(0) = (A'q)^ + (4.16 ± 0.50) x 10 



(2.10) 
(2.11) 



Ao" = (A^,)=^ - 6 ml/tcT G"(0) A^, - Sm^/t^T ^'(0) 

= (Aq)^ + 3 (4.16 ± 0.50) X 10"^ A'o + (2.72 ± 0.11) x 10"^ 



(2.12) 



Furthermore, thanks to Eq. ( p. 8]) , measuring /o(tcT) provides a significant constraint on 
JkI jvl f+ (0) limited only by the small theoretical uncertainty on Act- As we will discuss in 
Section ^.2.3| , this represents a powerful consistency check of present lattice QCD estimates 
of fx/fTT and /+(0). 

A similar dispersive parametrization for the vector form factor has been proposed in 



Ref. 1 11] with two subtractions performed at t = 0. This leads to: 



f+{t) = exp (A+ + H{t)) , where H{t) 



mtt 



(pis) 



°° ds 
tK. s^{s-t- ie) 



(2.13) 



In the elastic region, the phase of the vector form factor, (p{s), equals the / = 1/2, P-wave 
Kit scattering phase. 

Additional tests can be performed using the expression for the scalar form factor fo{t) 
at order / in ChPT M]: 



Ut) - 

where 
/+(0) 



/+(0) + A(t) + ^^^/^" }^ t + 4 (2Cr2 + CI,) [ml + ml)t - , (2.14) 



m 



K 



J TT 



Jn 



X'a = 8 



^0 



l + ^{^)- ^^{Cl^ + Cl,){ml-ml)' 

Jn 

ml {ml + ml) , ^ , 



n /+(o) 

,4 



m 



K 



K 



1 



1 



m2VA/+(0) /+(0) 



+ m. 



(2.15) 

" /+(0) 



-16 



+ < 



4A"(0) 



/+(0) 



Here A(t) is a function which receives contributions from order and p^, but like A(0) it 
is independent of the CJ, and the order p'* chiral constants L\ only appear at order p^. A{t) 
and A(0) have been evaluated in the physical region in Ref. ||l^ using for the L[ values a 
fit to experimental data. An analysis has been presented in ref. [|l^. However, the fit has 
to be reconsidered in light of the new experimental results as for instance considering the 
new Ki, analysis from NA48 and the updated value of fx/fn- 
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2.2.2 Analyticity and improved series expansion 

Armed only with the knowledge that the form factor is analytic outside the cut on the real 
axis, analyticity provides powerful constraints on the form factor shape without recourse 
to model assumptions. In particular, by an appropriate conformal mapping, the series 
expansion (2.5) necessarily "resums" into the form 

/(f) = i(ao + 01^ + 02^2 + ...), (2.16) 

where (p is an analytic function and 



,^t, to) = ^^2^.^^^ (2.17) 

is the new expansion parameter. In this "z expansion", the factor z{t,tQ) sums an infinite 
number of terms, transforming the original series, naively an expansion involving < 
0.3, into a series with a much smaller expansion parameter. For example, the choice 



to = tthi^ ~ V ~ ("^i: ~ '^'i)l^th) minimizes the maximum value of z occurring in the 
physical region, and for this choice |z(t,to)| ^ 0.047. 

The function (j) and the number to may be regarded as defining a "scheme" for the 
expansion. The expansion parameter z and coefficients at are then "scheme-dependent" 
quantities, with the scheme dependence dropping out in physical observables such as f{t). 
For the vector form factor, a convenient choice for (j) is 



4>F+{t,to,Q ) 



2, /I Z{t,0) ^z(t,-Q2)y/2 



327r -t V -Q - 1 

z{t,to)\'^^^ f z(t,t_)\'''^^'^ t+-t 



(2.18) 



This choice is motivated by arguments of unitarity, whereby the coefficients can be bounded 
by calculating an inclusive production rate in perturbation theory |2^. In fact, a much 
more stringent bound is obtained by isolating the exclusive Kir production rate in the 
vector channel from r decay data |2^. This enforces [pc| ] 

oo 9 



^^<170. (2.19) 

i.-n 



k=0 

With this choice of (f), and = 2 GeV^, a convenient choice for Iq is to = 0.39 {rriK — ^■n)'^ ■ 
This choice eliminates correlations in shape parameters ai/ag and a^ja^. 

The bound on the expansion coefficients can be used to bound errors on physical 
quantities describing the form factor shape, as discussed below in Sect. p.5| . A similar 
expansion can be used for the scalar form factor. Note that error estimates based on ( 2.ig| ) 



are conservative — no single coefficient is likely to saturate the bound. Also, this bound is a 
maximum taken over different schemes; more stringent bounds for particular schemes can 



be found in |2C]. 
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In addition to the direct applications in K^^ decays, it is important for other purposes 



to constrain the first few coefficients in (2.16), and check whether the series converges as 



expected. i^T^s decays provide a unique opportunity to do this. For example, the same 
parameterization can be used to constrain the form factor shape in lattice calculations 
of /(O), with the threshold tt/i adjusted to the appropriate value for the simulated quark 
masses. Measurements of in the kaon system can similarly be used to confirm scaling 
arguments that apply also in the charm and bottom systems pl|]. 

2.3 and K^2 decays beyond the SM 
2.3.1 The s ^ u effective Hamiltonian 

On general grounds, assuming only Lorentz invariance and neglecting effective operators 
of dimension higher than six, A5 = 1 charged-current transitions are described by 10 
independent operators: 

Ht'=' = -^Vus [clL{si^Lu){vrm + cli,{s^^Lu){u^^m) 

+ cidsrRu){yrLi) + cy,R{srRn){urm 

+ 4l{sLu){uM) + clji{sLu){i?m) 
+ cl^L{sRu){vLt) + cl^j^{sRu){vm) 

+clL{sa'"'Lu){ua'"'Le) + c^jiisa'"' Ru){ua'"' Ri)] + h.c. (2.20) 

where L = (1 — 75) and i? = (1 + 75). Defining this Hamiltonian at the weak scale, the SM 
case corresponds to cY^^{M^^) = 1 and all the other coefficients set to zero. The universal 
electromagnetic correction factor S'ew appearing in Eq. (|2.1[) describes the evolution of c^^^^ 
to hadronic scales: cYj^{M^)/cY^j^{M^) = 1 + (Sew - l)/2 f» . A similar expression can 
also be written for the Hamiltonian regulating u ^ d transitions. 

In the case oi K ^ iv decays only six independent combinations of these operators 
have a non- vanishing tree-level matrix element: 



A{K ^ ^£1/) = ( T^^y 



■;^cs(s'u)(i^^) + *]y:^cp(su)(z^75£) 



/If 2 ^ ' M2 '''5 



K 



(2.21) 



where 



cv = +(cL + + + cIr) , (2.22) 
CA = -{cIl + cIl - cIr - cIr) , (2.23) 
cs = +{cll + 4l + cIr + 4,R)Mw/mi , (2.24) 
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ICp 



crl 



clr 



CRR)Mw/mi , 



Similarly, in the K ^ case the independent structures are 



■^■^r7-^s{s-i-,u){v€) + ■^^kp{sj5u){uj5i) + h.c. 



where 



Mw 

ky 
kA 
ks 
kp 



Mw 



K 



V I V 
'^RL + cf _R 



cfij + 4R)Mw/mi . 



(2.25) 
(2.26) 



(2.27) 



(2.28) 
(2.29) 
(2.30) 
(2.31) 



On general grounds, new degrees of freedom weakly coupled at the scale Anp are ex- 



AS=1 



Focus- 



pected to generate corrections of ©(M^/A^p) to the Wilson coefficients of H 
ing on well- motivated new-physics frameworks, the following two scenarios are particularly 
interesting: 

• In two Higgs doublet models of type-II, such as the Higgs sector of the MSSM, 
sizable contributions are potentially generated by charged-Higgs exchange diagrams 

These are well described by the following set of initial 



(see e.g. Ref. ||, H, |§ 
conditions for s — > ti transitions. 



1 



and for u 



d transitions, 
1 



V,ud 



and 



and 



tan^/? 



minis 



"LR 



S,ud 

"lr 



(1 + eotan/3) 



tan^/? 



(1 + eo tan /?) 



(2.32) 



(2.33) 



H+ 



Here tan /3 is the ratio of the two Higgs vacuum expectation values and eg is a loop 
function whose detailed expression can be found in Ref. [^]. In presence of sizable 
sources of lepton-flavor symmetry breaking, a non-vanishing scalar-current contribu- 
tion to the lepton-flavor violating process K ev^- is also present ||2^. The latter 
can be parametrized by 

4^ tan^ /3 . (2.34) 



^LR 



m 



H+ 



In the Higgs- less model of Ref. , non-standard right-handed quark currents could 
become detectable. These are described by the following set of initial conditions for 
both — > s and u — > d transitions 



V,ud 



--{1 + 5) and 4i = e, , (2.35) 

"-LL - (1 + and c^^l'^ = Ens , (2.36) 
where Ex and d are free parameters of the model, can reach a few percents if the 
hierarchy of the right-handed mixing matrix is inverted. 
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2.3.2 Ke2 rates 

According to the Hamiltonian of Eq. ( p.27 ), the K12 rate of Eq. ( p.2| ) can be modified as 



^(-^£2(7)) _ Vus fxryiK fl- mym\^ ^ 



Vud 



f^rriT, V 1 - ml/ml 



X (1 + <5c 



(2.37) 



\kA - m\/{msMw)kp\'^ + \kv + m\/{msMw)ks\'^ 
^ \k''/ - ml/{mMw)kf\'^ + \k^'^ + ml/{mMw)kf\^ ' 

where m = rriu + rrid and kJj:'^ are defined for the u ^ d transition. In the MSSM scenario 



where 



tan^ [5 



H 



1 



1 + eo tan /? 



(2.38) 
(2.39) 



2.3.3 rates and kinematical distributions 



In the Ki^ case the non-standard operators of Eq. ( 2.21 ) could in principle modify the 
Daliz plot distribution. However, as we will show in the following, this effect turns out to 
be hardly detectable for most realistic new-physics scenarios. 

The hadronic form factors needed in the general case are the two FFs defined in 
Eq. ( |2.4D plus a tensor FF, whereas fo{t) allow us to parametrize also the scalar-current 
matrix element. More specifically, we have 



(vr- {k)\i-su)\K'{p)) = - 
{tt- {k)\{sa^''u)\K^{p)) = i 



m 



K 



m: 



{iris — ITT'^, 

pi^k" - p^'k^' 



rriK 



-Jo{t) , 
Bt [t) . 



(2.40) 
(2.41) 



The tensor form-factor was studied on the lattice 111, with the result Brit) ^ 1.2(1) /+ (0) /(l- 
0.3(l)i) at ^ ~ 2GeV in the MS scheme (an earlier order-of-magnitude estimate may be 
found in Ref.[|l4|). 

Choosing as independent kinematical variables 



'^PK ■ Ptt _ m\ + ml-t 



m 



K 



m 



y 



2PK 



mz 



K 



m 



K 



m 



K 



the double differential density can be written as (neglecting long-distance electromagnetic 
corrections) 



dr 



GlW, 



p\ "'tis I 



12^5 



m 



K 



CrSc 



dy dz 2567r^ 

+ As{y,z) (|5|2 + |P|2; 



Aiiy, z) (|F|2 + I - ^2(y, z) Re {VS* - AP* 



(2.42) 
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whereas 

Ai{y, z) =A{z + y- 1)(1 - y) + ri{Ay + 3z - 3) - 4r„ + r^(r^ - r^), 

A2{y,z) =2rt{?,-2y- z + rt-r^), A^iy, z) = re{l + - z - re). (2.43) 

Here S, P, V, and A are convenient combinations of hadronic form factors and short- 
distance Wilson coefficients: 

Vit,y) = Uit)cy-mj^'-I^ (2.44) 

A{t,y) = U{t)c^ + me^—-— 

Sit, y) = -ifiit) - U^t))"^^^^ cy - [ml + ml{2-z-2y)) JI^^eMI 

Pit,y) = if,^it)-Uit))<^c^-{mj+mli2-z-2y)) ^^^^^^^^ 



ruK 



where 



m = /„{<) ( 1 + r^^^^M- ') e.p (Sl^ 'A^) (2.45) 

V [ms-mu)Mw ) \ ly^w ms-muj 

h'(t) = m) ( 1 - t) . Mt) exp ( Zlilt±<^f" (2.46) 

V {ms-mu)Mw J V J^W nis-muj 

tcT = {ttik ~ ^n) '^^ have assumed cs,p/cv,a ^ 1- The SM case is recovered from 
Eq. (|2.44|) in the hmit cy = -ca = 1 and f^'^it) = fo{t). 



After integrating over y, differences to the SM rate of Eq. ( |2.1| ) can be summarized as 
it follows. Right-handed currents can only rescale the overall rate of Eq. (^]^), namely 

I 1 2 I 1 2 

r(i^,3(^)) ^ r(if,3(^)) X . (2.47) 

Scalar and pseudoscalar contributions can be easily encoded in Eq. ( |2.1|) by substituting 

In particular, these new effects are vanishing for i = 0, namely /o(0) in Eq. ( p.l| ) is free 
from them. The tensor coupling modify the phase space integral /|^(A-(-^o) of Eq. ( |2.1| ) by 



/^(A^,o) - /i(A-,,o) - ^^^^%^^^4(At,+,o) (2.49) 
In conclusion, the integrated rate including electromagnetic corrections can be written as 

(^2 ™5 / \ 2 

'^F^Kn..Q IT/ \2f /n\2 , xK , ^ 



r(i^^3(7)) = -^Ck Sew |Ks|^/+(0)^ (^1 + 6^ui2) + Cm j (2-50) 

IcyP + lcAp ( Re(c^4)-Re(c^^4) „\ 

X -iff- : tt; : rr; 1' 



11 



where 




i'k = z:t^ Idtxnt) (l + ^)(l-^) (2.51) 

flit) + 



dtXnt)^(l+^-^)(l-^X (2.52) 



' m2^/+(0)2 J " Amx V ' t J t 

BtH) X +6 ^t + 2mj)mj,Xit) ^'^'^ 

and A(t) = 1 - 2r2 + - 2t/mj^ - 2rlt/m\ + t^/m\. 

In most realistic new-physics scenarios the modification of the K^^ scalar form factor 
is well below the present experimental and theoretical errors. For instance, in the MSSM 
(or two-Higgs doublets) case f^{t) reads 

/o''(*)Imssm = foii) exp i-r^f'^^ (2.53) 

where is the parameter controlling the corrections to the rate of Eq. ( |2.38| ). For 
natural values of the free parameters (eo = 10~^, Mjj+ = 400 GeV and tan/3 = 40), such 
that = 0.2%, the corresponding modification of the Ki^ scalar form factor is 

^ ~ 1-0% or ^^^^^^^sM - 0-18% , (2.54) 
well below the level of present theoretical and experimental uncertainties. 
3. Data Analysis 

We perform fits to world data on the BRs and lifetimes for the Kl and K^, with the 
constraint that BRs add to unity. This is the correct way of using the new measurements. 
The fit procedure is described in Appendix 

3.1 Kl leading branching ratios and tl 

Numerous measurements of the principal BRs, or of various ratios of these BRs, have 
been published recently. For the purposes of evaluating 1 | /-|_ (0) , these data can be used 
in a PDG-like fit to the Kj^ BRs and lifetime, so all such measurements are interesting. A 
detailed description to the fit of the principal Kl BRs and tl is given in Appendix p| 



KTeV has measured five ratios of the six main Kl BRs |27]. The six channels in- 
volved account for more than 99.9% of the Kl width and KTeV combines the five mea- 
sured ratios to extract the six BRs. We use the five measured ratios in our analysis: 
BR(i^^3)/BR(ire3) = 0.6640(26), BR(7r+7r-7rO)/BR(iCe3) = 0.3078(18), BR(7r+^-)/BR(Ke3) 
0.004856(28), BR(37rO)/BR(iCe3) = 0.4782(55), and BR(27r°)/BR(37r°) = 0.004446(25). 
The errors on these measurements are correlated; this is taken into account in our fit. 
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I—' O "V O ^"1" tHY* 


Value 


o 






1.1 






1 1 
1 . 1 


-1— » J- i, V (J /I / 




1.2 


BR(7r+7r-7r°) 


0.1254(6) 


1.1 


BR(7r+7r-) 


1.997(7) X 10-3 


1.1 


BR(27rO) 


8.64(4) X 10"^ 


1.3 


BR(77) 


5.47(4) X 10-"^ 


1.1 


TL 


51.17(20) ns 


1.1 



Table 2: Results of fit to BRs and lifetime. 



NA48 has measured the ratio of the BR for K^^ decays to the sum of BRs for all 
decays to two tracks, giving BR(Ke3)/(l - BR(37r°)) = 0.4978(35) [||]. From a separate 
measurement of BR(i^L ^ 3tt^)/BR{Ks 2tt°), NA48 obtains BR(37r°)/rL = 3.795(58) 
/xs-i |2|. 

Using (j) — > KlKs decays in which the Ks decays to 7r"'"7r^, providing normalization, 
KLOE has directly measured the BRs for the four main decay channels |3C]. The 
errors on the KLOE BR values are dominated by the uncertainty on the Kl lifetime tl; 
since the dependence of the geometrical efficiency on is known, KLOE can solve for 
by imposing BR(^L — > x) = 1 (using previous averages for the minor BRs), thereby 
greatly reducing the uncertainties on the BR values obtained. Our fit makes use of the 
KLOE BR values before application of this constraint: BR^Kes) = 0.4049(21), BR(i^^3) = 
0.2726(16), BR(i^e3) = 0.2018(24), and BR{Ke3) = 0.1276(15). The dependence of these 
values on r/, and the correlations between the errors are taken into account. KLOE has 
also measured tl directly, by fitting the proper decay time distribution for Kl 3it^ 
events, for which the reconstruction efficiency is high and uniform over a fiducial volume 
of ~0.4Al. They obtain tl = 50.92(30) ns 

There are also two recent measurements of BR(7r+7r^)/BR(i^£3), in addition to the 
KTeV measurement of BR(7r^7r-)/BR(Ere3) discussed above. The KLOE collaboration ob- 
tains BR(7r+7r-)/BR(i^^3) = 7.275(68) x lO'^ |2|, while NA48 obtains BR(7r+7r-)/BR(A"e3) 
= 4.826(27) X 10^3 [^3|. All measurements are fully inclusive of inner bremsstrahlung. 
The KLOE measurement is fully inclusive of the direct-emission (DE) component, DE 
contributes negligibly to the KTeV measurement, and a residual DE contribution of 0.19% 
has been subtracted from the NA48 value to obtain the number quoted above. For consis- 
tency, in our fit, a DE contribution of 1.52(7)% is added to the KTeV and NA48 values. 
Our fit result for BR(7r^7r~) is then understood to be DE inclusive. 

In addition to the 14 recent measurements listed above, our fit for the seven largest 
Kl BRs and lifetime uses four of the remaining five inputs to the 2006 PDG fit and the 
constraint that the seven BRs add to unity. The results are given in Table ^. 

The evolution of the average values of the BRs for Kf^i^ decays and for the important 
normalization channels is shown in Fig. p|. 
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Figure 1: Evolution of average values for main BRs. 



Our fit gives xV^df = 20.2/11 (4.3%), while the 2006 PDG fit gives xV^df = 14.8/10 
(14.0%). The differences between the output values from our fit and the 2006 PDG fit are 
minor. The poorer value of x^/ndf for our fit can be traced to contrast between the KLOE 
value for BR(37r'') and the other inputs involving BR(37r'^) and BR(7r''7r'') — in particular, 
the PDG ETAFIT value for BR(7r°7r°/7r+7r"). The treatment of the correlated KTeV and 
KLOE measurements in the 2006 PDG fit gives rise to large scale factors for BR(i^e3) and 
BR(37r'^); in our fit, the scale factors are more uniform. As a result, our value for BR(^Ce3) 
has a significantly smaller uncertainty than does the 2006 PDG value. 



3.2 Ks leading branching ratios and rg 



KLOE has published a measurement of BR(Er5 — > irev) that is precise enough to 
contribute meaningfully to the evaluation of |V^s|/+(0). The quantity directly measured 
is BR(7rez/)/BR(7r+7r^). Together with the published KLOE value BR(7r+7r-)/BR(7r°7r°) 
= 2.2459(54), the constraint that the Ks BRs must add to unity, and the assumption of 
universal lepton couplings, this completely determines the Ks BRs for vr+vr", vr'^vr'^, K^z, 
and ET^s decays |5|. In particular, BR(if5 ireu) = 7.046(91) x 10"^. 

NA48 has recently measured the ratio T{Ks ■Kev)/T{KL vrei^) = 0.993(26) (22) 
1 36]. The best way to include this measurement in our analysis would be via a combined 



fit to Ks and K^ branching ratio and lifetime measurements. Indeed, such a fit would be 
useful in properly accounting for correlations between Ks and Kl modes introduced with 
the preliminary NA48 measurement of T{Kl — > 37r^), and more importantly, via the PDG 
ETAFIT result, which we use in the fit to Kl branching ratios. At the moment, however, 
we fit Ks and Kl data separately. NA48 quotes BYi{Ks Treu) = 7.046(180) (160) x 10"^; 
averaging this with the KLOE result gives BIl{Ks — > vrez/) = 7.05(8) x 10^^, improving 
the accuracy on this BR by about 10%. 

For we use 0.8958 x 10^^'^ s, where this is the non-CPT constrained fit value from 
the PDG, and is dominated by the 2002 NA48 and 2003 KTeV measurements. 



3.3 K leading branching ratios and r 

There are several new results providing information on K^ rates. These results are mostly 
preliminary and have not been included in previous averages. 
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1 1 
1 . 1 


BRfTTTr^'l 

J— 'J- L I /I / 1 J 


20 64('8)% 


1.1 


BR(7r7r7r) 


5.595(31)% 


1.0 


BR(i^e3) 


5.078(26)% 


1.2 




3.365(27)% 


1.7 


BR(7r7r°7rO) 


1.750(26)% 


1.1 


r± 


12.384(19) ns 


1.7 



Table 3: Results of fit to BRs and lifetime. 



The NA48/2 collaboration has recently published measurements of the three ratios 



BR(i^'e3/7r7r*'), BR^K^^/mr^), and BR{K^3/ K^^) [37|. These measurements are not in- 
dependent; in our fit, we use the values BK{K(,3/m:^) = 0.2470(10) and BR(il'^3/7r7r'^) = 
0.1637(7) and take their correlation into account. ISTRA+ has also updated its preliminary 
value for BR(i^e3/7r7r°). They now quote BR(i^e3/7r7r°) = 0.2449(16) |8|]. 

KLOE has measured the absolute BRs for the K^s and K^^^ decays [ p9| . In </> — > K~^K~ 
events, decays into fiv or tttt'^ are used to tag a K~ beam, and vice versa. KLOE per- 
forms four separate measurements for each Kf^^ BR, corresponding to the different combina- 
tions of kaon charge and tagging decay. The final averages are BR(^Ce3) = 4.965(53)% and 
BR(i^^^3) = 3.233(39)%. Very recently KLOE has also measured the absolute branching 
ratio for the tttt^ decay with 0.5% accuracy. The KLOE preliminary result, is BR(7r7r'') = 
0.20658(112) ll]. 

Our fit takes into account the correlation between these values, as well as their depen- 
dence on the lifetime. The world average value for t± is nominally quite precise; the 
2006 PDG quotes t± = 12.385(25) ns. However, the error is scaled by 2.1; the confidence 
level for the average is 0.17%. It is important to confirm the value of t±. The two new 
measurements from KLOE, t± = 12.367(44) (65) ns and t± = 12.391(49)(25) nsQ with 
correlation 34%, agree with the PDG average. 

Our fit for the six largest BRs and lifetime makes use of the results cited above, 
plus the data used in the 2006 PDG fit, except for the Chiang '72 measurements for a total 
of 26 measurements. The six BRs are constrained to add to unity. The results are shown 
m Table |. 

The fit quality is poor, with x^/ndf = 42/20 (0.31%). However, when the five older 
measurements of t± are replaced by their PDG average with scaled error, x^/ndf improves 
to 24.3/16 (8.4%), with no significant changes in the results. 

Both the significant evolution of the average values of the Ki^ BRs and the effect of 
the correlations with BR(7r7r'^) are evident in Fig. |2[ 

3.4 Measurement of BR(E:e2)/BR(K^.2) 

Experimental knowledge of Kf>2/ K^2 has been poor so far. The current world average of 
Rk = BR(iv:e2)/BR(i^^2) = (2.45 ±0.11) x 10"^ dates back to three experiments of the 
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Figure 2: Evolution of average values for main BRs. 



1970s |j43| and has a precision of about 5%. Three new prehminary measurements vi^ere 
reported by NA48/2 and KLOE (see Tab. |): A prehminary resuh of NA48/2, based on 
about 4000 Ke2 events from the 2003 data set, was presented in 2005 [^4|. Another pre- 
liminary result, based on also about 4000 events, recorded in a minimum bias run period 
in 2004, was shown at KAONOT fl^ . Both results have independent statistics and are also 
independent in the systematic uncertainties, as the systematics are either of statistical na- 
ture (as e.g. trigger efficiencies) or determined in an independent way. Another preliminary 
result, based on about 8000 i^e2 events, was presented at KAON07 by the KLOE collabo- 
ration Both, the KLOE and the NA48/2 measurements are inclusive with respect to 
final state radiation contribution due to bremsstrahlung. The small contribution of Ki2-y 
events from direct photon emission from the decay vertex was subtracted by each of the 
experiments. Combining these new results with the current PDG value yields a current 
world average of 

Rk = (2.457 ± 0.032) x 10"^ (3.1) 

in very good agreement with the SM expectation and, with a relative error of 1.3%, a factor 
three more precise than the previous world average. 

3.5 Measurements of K^^ slopes 

3.5.1 Vector form factor slopes from K^^ 

For Ke3 decays, recent measurements of the quadratic slope parameters of the vector form 
factor (A^, A'|) are available from KTeV |4|, KLOE ||, ISTRA+ [||, and NA48 g. 
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Rk [10-5] 


PDG 2006 ||] 


2.45 ±0.11 


NA48/2 prel. ('03) || 


2.416 ± 0.043 ± 0.024 


NA48/2 prel. ('04) || 


2.455 ± 0.045 ± 0.041 


KLOE prel. |§ 


2.55 ±0.05 ±0.05 


SM prediction 


2.477 ±0.001 



Table 4: Results and prediction for Rk = BR{Ke2)/^R{K^2)- 



We show the results of a fit to the and K~ data in the first column of Table ^, and 
to only the Kl data in the second column. With correlations correctly taken into account, 
both fits give good values of x'^/ndf. The significance of the quadratic term is 4.2(7 from 
the fit to all data, and 3.5a from the fit to Kl data only. 





Kl and K data 
4 measurements 
xVndf = 5.3/6 (51%) 


Kl data only 
3 measurements 
xVndf = 4.7/4 (32%) 


A'+ X 10^ 


25.2 ±0.9 


24.9 ± 1.1 


A'; X 10^ 


1.6 ±0.4 


1.6 ±0.5 


P(AV,A±") 


-0.94 


-0.95 




0.15465(24) 


0.15456(31) 




0.15901(24) 


0.15891(32) 



Table 5: Average of quadratic fit results for K^^ slopes. 



Including or excluding the K~ slopes has little impact on the values of A^ and \"^; 
in particular, the values of the phase-space integrals change by just 0.07%. The errors on 
the phase-space integrals are significantly smaller when the K~ data are included in the 
average. 

KLOE, KTeV, and NA48 also quote the values shown in Table ^ for My from pole fits 
to Kl e3 data. The average value of My from all three experiments is My = 875 ± 5 MeV 
with x^/iidf = 1.8/2. The three values are quite compatible with each other and reasonably 
close to the known value of the K=t*(892) mass (891.66 ± 0.26 MeV). The values for A',_ 
and \'\_ from expansion of the pole parametrization are qualitatively in agreement with the 
average of the quadratic fit results. More importantly, for the evaluation of the phase-space 
integrals, using the average of quadratic or pole fit results gives values of I{K^'^) that differ 
by just 0.03%. 

An attempt to estimate the theoretical uncertainties associated to form factor parame- 
terization has been pursued by KTeV, analyzing Kf.'i decays with the z-expansion 2.2. I for 



the f+{t) form factor g^]. The results are ai/oo = 1.023 ± 0.040 and 02/00 = 0.75 ± 2.16. 
The second order term is consistent with zero and the higher orders are bounded by the 
theory: X^^g '^fc/'^o — ^'''O- Using these results the phase space integral is calculated to be 
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Experiment 


Mv (MeV) 


(My) = 875 ± 5 MeV 


KLOE 


870 ± 6 ± 7 


xVndf = 1.8/2 


KTeV 


881.03 ±7.11 


A'^_ X 10^ = 25.42(31) 


NA48 


859 ± 18 


A'l = 2 X 
I{K%) = 0.15470(19) 



Table 6: Pole fit results for K^^ slopes. 



/(isrOg) = 0. 15392 ± 0.00048oxp ± 0.00006th. The first error corresponds to the KTeV exper- 
imental uncertainty and the second error is due to possible effects from higher order terms 
in the 2;-expansion. Compared to the global average using the quadratic parameterization 
(Table 5), the KTeV measurement using the z-expansion deviates by about 1.5(Texp. This 
result is less precise statistically, but it is more conservative as far as the estimate of the 
theoretical uncertainty is concerned. 

3.5.2 Scalar and Vector form factor slopes from ^^^3 

For i^^3 decays, recent measurements of the slope parameters (A'|_, A", Aq) are available 



from KTeV [|§, KLOE [|2[, ISTRA+ |5|1, and NA48 Q. These data are summarized in 
Appendix |D[ 

We have studied the statistical sensitivity of the form-factor slope measurements using 
Monte Carlo techniques, see Appendix]^.. The conclusions of this study are a) that neglect- 
ing a quadratic term in the parameterization of the scalar form factor when fitting results 
leads to a shift of the value of the linear term by about 3.5 times the value of the quadratic 
term; and b) that because of correlations, it is impossible to measure the quadratic slope 
parameter from quadratic fits to the data at any plausible level of statistics. The use of 
the linear representation of the scalar form factor is thus inherently unsatisfactory. 

Figure ^ shows the 1-a contours from all the experimental results (-f^es and K^-^). It is 
immediately clear from the figure that the new NA48 results are difficult to accommodate 
^ . Performing the combination with and without the NA48 results for the -R'^a form- factor 
slopes included we obtain fit probability values of 1 x 10~^ and 22.3% respectively(see 
Appendix ^ for a detailed comparison). The results of the combination are listed in 
Table |^ 

The value of x^/ndf for all measurements is terrible; we quote the results with scaled 
errors. This leads to errors on the phase-space integrals that are ~60% larger after inclusion 
of the new JC^s NA48 data. 

We have checked to see if the NA48 data might show good consistency with the 
results of some other experiment in a less inclusive average. Fitting to only the results 
from KTeV, NA48, and ISTRA+ gives xV^df = 28/6 (0.01%). Fitting to only the KL^,■i 
results from KTeV, NA48 gives xV^df = 12/3 (0.89%). The consistency of the NA48 data 
with these other measurements appears to be poor in any case. 

The evaluations of the phase-space integrals for all four modes are listed in each case. 
Correlations are fully accounted for, both in the fits and in the evaluation of the integrals. 

^It lies out of correlation directions in the [A^, A" , Ao] space 
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KLOE KTeV ISTRA+ NA48 




Figure 3: l-a contours for A^, A" , Aq determinations from ISTRA+(pink ellipse), KLOE(blue 
ellipse), KTeV(red ellipse), NA48(green ellipse), and world average with(filled yellow ellipse) and 
without (filled cyan ellipse) the NA48 result. 



Kl and K Kl only 

Measurements 16 11 

xVndf 54/13 (7 x 10"'^) 33/8 (8 x 10"^) 

AV X 10^ 24.9 ±1.1 (5 = 1.4) 24.0 ±1.5 (5 = 1.5) 

A'l X 10^ 1.6 ± 0.5 (S = 1.3) 2.0 ± 0.6 (S = 1.6) 

Ao X 103 13.4 ±1.2 (5 = 1.9) 11.7 ± 1.2 (5 = 1.7) 

p(A'+,A'|) -0.94 -0.97 

p(A'+,Ao) ±0.33 ±0.72 

p(A'|,Ao) -0.44 -0.70 

I{K%) 0.15457(29) 0.1544(4) 

liKf^) 0.15892(30) 0.1587(4) 

I(i^°3) 0.10212(31) 0.1016(4) 

/(ivrj) 0.10507(32) 0.1046(4) 

p{Ie3,I^3) ±0.63 ±0.89 



Table 7: Averages of quadratic fit results for K^s and i^T^a slopes. 

The correlation matrices for the integrals are of the form 

±1 ±1 p p 
±1 ±1 p p 

p p ±1 ±1 

p p ±1 ±1 

where the order of the rows and columns is K^^, K^, K'^^, K^^, and p = p{Ie3, I^s) as 
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0.28 
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0.19 


0.15 


0.09 


Kl ^ 


TTflU 


0.2168(7) 


0.31 


0.10 


0.18 


0.15 


0.15 


Ks^ 




0.2154(13) 


0.67 


0.65 


0.03 


0.15 


0.09 


K^- 


> nev 


0.2173(8) 


0.39 


0.26 


0.09 


0.26 


0.09 


K^- 




0.2176(11) 


0.51 


0.40 


0.09 


0.26 


0.15 


averag 




0.2166(5) 













Table 8: Summary of \Vus\ x /+(0) determination from all channels. 



listed in the table. 

Adding the data to the fit does not cause drastic changes to the values of the 
phase-space integrals for the Kgs modes: the values for I{K^^) and I{K^) in Table |^ are 
qualitatively in agreement with those in Table ^. As in the case of the fits to the i^es data 
only, the significance of the quadratic term in the vector form factor is strong (3.6(7 from 
the fit to all data). 

4. Physics Results 

4.1 Determination of |Ks| x /+(0) and |T4s|/|Kd| x Jk/J^ 

This section describes the results that are independent on the theoretical parameters /+(0) 
and fK/fn- 

4.1.1 Determination of x /+(0) 



The value of \Vus\ x /+(0) has been determined from (2J) using the world average values 
reported in section |^ for lifetimes, branching ratios and phase space integrals, and the 
radiative and SU (2) breaking corrections discussed in section ||. 

The results are given in Table |8|, and are shown in Fig. |^ for Kl ireu, Kl tt/uz/, 
Ks — > TTeu, K — > Treu, K^ — > TTiUi/, and for the combination. The average, 

|V;.| X /+(0) = 0.21664(48), (4.1) 

has an uncertainty of about of 0.2%. The results from the five modes are in good agreement, 
the fit probability is 58%. In particular, comparing the values of \Vus\ x /+(0) obtained 
from K^^ and K^ we obtain a value of the SU(2) breaking correction 

^sui2U.. = 2-9(4)% 

in agreement with the CHPT calculation reported in Table ||: <5_^[/(2) — 2.36(22)%.^ 



^ The value of S§^^2) has a direct correspondence to the ratio of light quark masses. Recent analyzes 
on the so-called violations of Dashen's theorem in the Kaon electromagnetic mass splitting point to Sgui^2) 
values of about 3% |56|. 
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Figure 4: Display of \Vus\ x /+(0) for all channels. 



4.1.2 Determination of |Ks|/|Kd| x /^-//^ 

An independent determination of \Vus\ is obtained from decays. The most impor- 
tant mode is fJ'^'^, which has been recently updated by KLOE reaching a relative 
uncertainty of about 0.3%. As shown in Eq. (|2.2|), hadronic uncertainties are minimized 
considering the ratio T{K~^ fi'^ u) /T {tt'^ f^^^)- 

Using the world average values of BR(K=t fs^y) and of given in Section ^ and 
the value of r(7r=^ fi^u) = 38.408(7) j^s'^ from |4|] we obtain: 

\Vus\/\Vud\ X fx/f. = 0.2760 ± 0.0006 . (4.2) 



4.2 The parameters /+(0) and /k/ fir 

The main obstacle in transforming these highly precise determinations of \Vus\ x /+(0) and 
I I /I Kid I X JkI f-K into a determination of |Kis| at the per- mil level are the theoretical 
uncertainties on the hadronic parameters /+(0) and Jk/ fw 



4.2.1 Theoretical estimates of /+(0) 

By construction, /+(0) is defined in the absence of isospin-breaking effects of both elec- 



tromagnetic and quark-mass origin. More explicitly, as discussed in Section 2.1, /+(0) is 
defined by the vr"*" matrix element of the vector current in the limit = rrid and 

ctem — ^ 0, keeping kaon and pion masses to their physical values. 

This hadronic quantity cannot be computed in perturbative QCD, but it is highly 
constrained by SU{3) and chiral symmetry. In the chiral limit and, more generally, in 
the SU{3) limit (m^ = rrid = nis) the conservation of the vector current (CVC) implies 
/-(_(0) = 1. Expanding around the chiral limit in powers of light quark masses we can write 

/+(0) = l + /2 + /4 + ... (4.3) 
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Figure 5: Present determinations of /+(0) = ff^ 
analytical approaches p^, B§, M, M . 



(0) from lattice QCD and analytical or semi- 



where /2 and /4 are the NLO and NNLO corrections in ChPT. The Ademollo-Gatto 
theorem implies that [/+(0) — 1] is at least of second order in the breaking of SU{3) or in 
the expansion in powers of — m, where m = {rriu + md)/2. This in turn implies that /2 
is free from the uncertainties of the 0{p'^) counterterms in ChPT, and it can be computed 
with high accuracy: /2 = —0.023 jSTjl . 

The difficulties in estimating /+(0) begin with /4 or at 0{p^) in the chiral expansion. 
At this order we can write 

h = ^{^i) + h\'-'{^i), (4.4) 

where A(;u) is the loop contribution, which has been computed in Ref. and /4p°'^(^) 
is the 0{p^) local contribution, whose knowledge cannot be simply deduced from other 
processes. Several analytical approaches to determine /4 have been attempted over the 
years [^], essentially confirming the original estimate by Leutwyler and Roos [57| (see 
Fig. |5|). The benefit of these new results, obtained using more sophisticated techniques, 
lies in the fact that a better control over the systematic uncertainties of the calculation 
has been obtained. However, the size of the error is still around or above 1%, which is not 
comparable to the 0.2% accuracy which has been reached for \Vus\ x /+(0). 

Recent progress in lattice QCD gives us more optimism in the reduction of the error on 
/+(0) below the 1% level [^, |6^. Most of the currently available lattice QCD results 
have been obtained with relatively heavy pions and the chiral extrapolation represents the 
dominant source of uncertainty. As shown in Figure ^, there is a general trend of lattice 
QCD results to be slightly lower than analytical approaches. An important step in the 



22 



TTTTTTTTTTTT-p-n-r 



N=0 



N^2 



N=2^1 



CP-PACS 

CP -P ACS 
JLQCD 
RBC 

QCDSF-07 
ETMC-07 

MILC-04 
MILC-07 
NPLQCD-07 

PACS-CS-07* 
RBC-UKQCD-07<j 
HPQCD-UKQCD-Cj7 



♦ 



a=0.05fm->0, m >500 MeV 



1.192(30) 


Wilson: a=0.l 


1fm->C 


), m > 


1.148," 


Wilson: 


m 


„>550 


MeV 


1.175(11) 






m 


^>550 


MeV 


1.219(26) 


Wilson: 


m 


^>300 


MeV 


1.227(26) 


TMQCD 


rr 


>300 

It 


MeV 



1.210(14) Stag. 



1.197 



13 Stag. 

24 



a=0.06 fm -> ;;;;;; 
m >280 MeV, Lm >5.0 
l-218ji'"' DWF/Stag. 
1.219(22) Wilson : m>2\ MeV, Lm^>2.9 
1.205(18) 



1.189(7) Stag. 



I I I I I I I I I I I I I I I I I 



v.v^ v^-* \^ 

Figure 6: Summary of JkI !tx estimates |6^, |6^. All values are from Lattice QCD. In 

recent studies, sea quarks are getting light and data are matched to ChPT fits to determine the 
Low-energy-Constants (LEG). 



reduction of the error associated to the chiral extrapolation has been recently made by the 
UKQCD-RBC cohaboration |8|. Their prehminary result /+(0) = 0.964(5) is obtained 
from the unquenched study with Np = 2 + 1 flavors, with an action that has good chiral 
properties on the lattice even at finite lattice spacing (domain- wall quarks). They also 
reached pions masses (> 330 MeV) much lighter than that used in previous studies of 
/+(0). The overall error is estimated to be 0.5%, which is very encouraging. Moreover, 
they observe for /+(0) a mass dependence similar to the one of f^- That is something new 
with respect to previous lattice studies (this is likely due to the fact that they work with 
lighter pions). To assess the chiral uncertainty of /4, polynomial fits (linear and quadratic) 



well reproduce the data. However, it would be interesting |61] to have the expression of 
A(/i) in Eq. (|4.4D in terms of the quark masses so to directly estimate /4|'°'^(/i)- Moreover, 
it should also be stressed that the present study is performed at a single value of the lattice 
spacing (a = 0.12 fm) and in a relatively small extension of the fifth dimension of the 
lattice.^ 

In the following phenomenological analysis we will use this result as the present best 
estimate of /+(0), although some reservation remains. 



^ Even though m^L > 4.5, simulations with a larger fifth dimension, Ls would help too because the 
mass of their lightest quark (= 0.005 in lattice units) is very close to the residual mass parameter 0.003, 
also in lattice units). This may entail some uncontrolled systematics, in particular for Jk/ .f-n 
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4.2.2 Theoretical estimates of Jk/ fir 

In contrast to the semileptonic vector form factor, the pseudoscalar decay constants are 
not protected by the Ademollo-Gatto theorem and receive corrections linear in the quark 
masses. Expanding JkI J-k in power of quark masses, in analogy to /+(0), 

fKlh = 1 + r2 + . . . (4.5) 

one finds that the 0(jf) contribution r2 is already affected by local contributions and 
cannot be unambiguously predicted in ChPT. As a result, in the determination of fK/f-w 
lattice QCD []64| - ||67|] has essentially no competition from purely analytical approaches. 
The status of the lattice results for fx/ f-K is summarized in Fig. |^. As can be seen, 
the present overall accuracy is about 1%. The novelty are the new lattice results with 
Np = 2 + 1 dynamical quarks and pions as light as 280 MeV |^, obtained by using 
the so-called staggered quarks."^ The analyzes of [|4|, cover a broad range of lattice 
spacings (i.e. a=0.06 and 0.15 fm) and is performed on sufficiently large physical volumes 
(uIt^L > 5.0). It should be stressed, however, that the sensitivity of fx/ fw to lighter pions 
is larger than in the computation of /+(0) and that chiral extrapolations are far more 
demanding in this case.^ Notice also that at Lattice 2007 preliminary studies with Np = 
2 + 1 clover quarks and pion masses > 200 MeV have been presented from either PACS-CS 
Collaboration ||7^ and ref. |7^. With respect to the results obtained with staggered quarks, 
the PACS-CS value of fx/ fw in fig- ^ is restricted to a single lattice spacing (a = 0.09 fm) 
and relatively small physical volume {m-jrL > 2.9). For ref. [^, the final analysis is to 
be completed. In the following analysis we will use as reference value the MILC-HPQCD 
result /k/ f-K = 1.189(7) [^], although some reservation about staggered fermions remains. 

4.2.3 A test of lattice calculation: the Callan-Treiman relation 

As described in Sect. the Callan-Treiman relation fixes the value of scalar form factor 
at t = — (the so-called Callan-Treiman point) to the ratio (/i<-//7r)//+(0). The 
dispersive parametrization for the scalar form factor proposed in |jl^ and discussed in 



* Staggered fermions come in four tastes on the lattice. In the continuum limit the extra degrees of 
freedom decouple from physical predictions. But, at finite lattice spacing, where the data are produced, 
the taste symmetry is violated and the extra degrees of freedom are removed by hand, namely by taking 
the fourth root of the staggered quark determinant. Theoretically, this procedure has been only confirmed 
in perturbation theory and is currently a subject of controversies within the lattice QCD community [^sj. 
Since the staggered dynamical quarks are computationally cheap, they have been largely used. Thanks 
to recent progress in algorithm building |^^, safer but still computationally competitive alternatives are 
becoming available. 

^ In some details, efi'ects of chiral logs are not clearly disentangled and analytic terms (NNLO or NNNLO) 
are still needed in order to extrapolate from the simulated sea quark masses (such as > 280 MeV) to 
the physical point. For example, the two studies of ref. and of ref. jHsj with staggered quarks share the 
same configurations, but they differ in how to extrapolate to the physical masses. Then, the central values 
of /x/Ztt between the two analyzes (namely, /k/Ztt ~ 1.197i3 s^nd fK/f-w = 1.189(7) from ref. j64j and 
ref. respectively) differ for la. However, taking into account the complete uncertainty of fx/ f-n in [ p^ , 
we have fK/ fw = 1.194(10) of j6^] for a symmetric error and the values of ref. and ref. jHsj look now 
in good agreement. The highly improved staggered fermions (HISQ) used in ^] for the valence quarks are 
designed to reduce the taste violation effects, which also should reduce the overall systematic uncertainty. 
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Figure 7: Values for f+ (0) determined from the scalar form factor slope using tlie Callan-Treiman 
relation and fK/fn = 1.189(7). The UKQCD/RBC result /+(0) = 0.964(5) is also shown. 



Sect. allows to transform the available measurements of the scalar form factor into a 
precise information on (/x//7r)//+(0), completely independent of the lattice estimates. 

Very recently KLOE []52| and NA48 ||5^ have presented results on the scalar FF slope 
using the dispersive parametrization. In these analyzes a dispersive parametrization is used 
for both the scalar and the vector form factors. A similar analysis has started for the KTeV 
data. We report these preliminary results for the first time. The ISTRA+ measurement of 
the scalar form factor slope performed using the first order Taylor expansion parametriza- 
tion can be translated in the dispersive parametrization as described in Appendix ^. The 
results are given in Table ^ for all the four experiments in the case of the pole parametriza- 
tion for the vector form factor. The original KLOE and NA48 results are also shown for 
comparison as well as the preliminary result of KTeV obtained from the K^^ data analysis. 
Moreover, a combined K^^ and K^^ data analysis is also in progress and the preliminary 
result is: log(C) = 0.191 it 0.012. The preliminary KTeV results are obtained using the 



original MC and data from Ref. |48] 



Experiment 


log(C) direct 


log(C)t 


KTeV 


0.195(14)* 


0.203(15) 


KLOE 


0.207(24) 


0.207(23) 


NA48 


0.144(14) 


0.144(13) 


ISTRA+ 




0.226(13) 



^ Estimated from Aq published. *Prcliminary results. 

Table 9: Experimental results for log(C). 

Figure ^ shows the values for /+(0) determined from the scalar form factor slope mea- 
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surements obtained using the Callan-Treiman relation and fK/f-w = 1-189(7). The value of 
/+(0) = 0.964(5) from UKQCD/RBC is also shown. As already noticed in Section |3|, the 
NA48 result is difficult to accommodate. Here one can see that this results is also not con- 
sistent with the theoretical estimates of /+(0). In particular, it violates the Fubini-Furlan 
bound /+(0) < 1 For this reason, the NA48 result will be excluded when using the 
Callan-Treiman constraint. 

The average of the experimental results on the FFs with the pole parametrization for 
the vector case and the dispersive parametrization for the scalar FF give: 

Ai = 0.0256 ± 0.0002 , 

(4.6) 

\l = 0.0149 ± 0.0007 , 

with correlation coefficient —0.32. The above results are then combined with the lattice 
determinations of fK/f-K = 1.189(7) and /+(0) = 0.964(5) using the constraint given by 
the Callan-Treiman relation. The results of the combination are given in Table |l^, where 
log C = \l tcT/ml + 0.0398 ± 0.0041. 



K 




/+(0) 


fK/f-K 


0.02563(19) 


0.0146(5) 


0.96(4) 


1.192(6) 


correlation matrix 


1. 


-0.23 


0.12 


-0.14 




1. 


-0.51 


0.61 






1. 


0.30 








1. 



Table 10: Results from the form factor fit. 



The fit probability is 39%, confirming the agreement between experimental measure- 
ments and lattice determination. The accuracy of fK/f-w is also slightly improved, and this 
effect can be better seen in the ratio f^{0)/{fK/f-K), directly related to the Callan-Treiman 
constraint. 

As previously discussed, new physics contributions to the scalar form factor (reab- 
sorbed into the value of logC) are generated only by scalar operators. Hence in the case 
of right-handed currents logC coincides with the SM value. This imply we can use the 
Callan-Treiman improved f+{0)/{fK/f-K) in constraining right-handed currents. On the 
other hand, this is not possible in the MSSM scenario, where scalar operators are present. 
Here the measured value of logC, following from ( 2.53| ), is 



logC^SSM ^ ^ i^g^SM _ (4^7) 



with the given in ( 2.39| ). By construction, the quantity logC^^ depends only on QCD 



dynamics and must satisfies the Callan-Treiman relation (^^. The theoretical calculation 
of /o(^ct) can thus be used to constrain scalar densities. At present, the theoretical 
knowledge of logC^^ is obtained from Eq. ( p.^ ) and is limited by our knowledge of Act, 
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Figure 8: Results of fits to |K.|, and 

reported in (^), and by the lattice QCD results on {fx/ fir)/ f+{0)- Using this information 
we obtain the constraint 

-0.007 ± 0.012 . (4.8) 



To improve this result it would be particularly useful a direct computation of (/x//7r)//+(0) 
on the lattice (i.e. from the the same set of simulations). Given the advanced status of 
staggered results on fx/ f-w, it would be interesting to see the effect of a corresponding 
analysis /+(0) (which at present is still very preliminary []59| ). 



4.3 Test of Cabibbo Universality or CKM unitarity 

To determine \ Vus\ and \ Vud\ we use the value |T4s| x /+(0) = 0.2166(5) reported in Table |8|, 
the result \Vus\/\Vud\fK / f-n = 0.2760(6) discussed in Sect. |4.1.2| , /+(0) = 0.964(5), and 
fK/f-K = 1.189(7). From the above we find: 



\Vus\ 
\Vus\/\Vud\ 



0.2246 ± 0.0012 
0.2321 ±0.0015 



[Ki^ only] , 
[Ki>2 only] . 



(4.9) 
(4.10) 



These determinations can be used in a fit together with the the recent evaluation of \Vud\ 
from 0+ 0+ nuclear beta decays: |Kd|=0. 97418 ± 0.00026 p|. The global fit gives 



\Vud\ = 0.97417(26) |K«| = 0.2253(9) [i^a,^2 + 0+ ^ 0+] 



(4.11) 



with x^/ndf = 0.65/1 (42%). This result does not make use of CKM unitarity. If the 
unitarity constraint is included, the fit gives 



\Vus\ = sin 9c = X = 0.2255(7) [with unitarity] 



(4.12) 
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s, d 

Figure 9: Z' and Higgs exchange 




and xVndf = 0.80/2 (67%). Both results are illustrated in Fig. |. 

As described in the introduction, the test of CKM unitarity can be also interpreted as 
a test of universality of the lepton and quark gauge couplings. Using the results of the fit 
(without imposing unitarity) we obtain: 

GcKM = [\Vud\^ + \Vus\^ + \Vub?]^'^ = (1.1662 ± 0.0004) x 10"^ GeV'^ , (4.13) 
in perfect agreement with the value obtained from the measurement of the muon lifetime: 



= (1.166371 ± 0.000007) x 10"^ GeV 



(4.14) 



The current accuracy of the lepton-quark universality sets important constraints on model 
building beyond the SM. For example, the presence of a Z' (see Fig. ^, left) would affect 
the relation between Gckm and G^ in the following way, 



Gil G: 



CKM 



l-Qm7QeL{Qt,L-Q 



dh) 2 1 2 T ' (^-^^^ 

where Qfi are the generic charges of the Z' to left-handed leptons (in units of the SM 
SU{2)l charge). In case of a Z' from SO{W) grand unification theories (QeL = Q^iL = 
—^QdL = 1) we obtain mz' > 700 GeV at 95% CL, to be compared with the mz' > 720 GeV 



bound set through the direct collider searches |43]. In a similar way, the unitarity constraint 



also provides useful bounds in various super symmetry-breaking scenarios |74]. 



4.3.1 Bounds on helicity-suppressed amplitudes 

A particularly interesting test is the comparison of the \Vus\ value extracted from the 
helicity-suppressed Ki2 decays with respect to the value extracted from the helicity-allowed 
Kg^ modes. To reduce theoretical uncertainties from fx and electromagnetic corrections 
in K£2, we exploit the ratio Br[Ki2) / Br[-Kf2) and we study the quantity 



Vus{Kn) Kd(0+ ^ 0+) 



(4.16) 



Vus{Kl^) Vud{l^t2) 

Within the SM, i?;23 = 1, while deviation from 1 can be induced by non- vanishing scalar- 
or right-handed currents. Notice that in i?/23 the hadronic uncertainties enter through 
(/x//.)//+(0). 



Following the notation of Section 2.5, effects of scalar currents due to a charged Higgs 
(Fig. P right) give 



123 



m 



K+ 



^ rrid 



tan^/3 



1 -|- eo tan (3 



(4.17) 
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Figure 10: Excluded region in the cliarged Higgs mass-tan /3 plane. The region excluded by 
B ^ Tiy is also indicated. 



whereas for right-handed currents we have 

Rl23 = 1 - 2 (e, - Ens) . (4.18) 

In the case of scalar densities (MSSM), the unitarity relation between \Vud\ extracted 
from 0^ — > 0^ nuclear beta decays and \Vus\ extracted from Kg^ remains valid as soon as 
form factors are experimentally determined. This constrain together with the experimental 
information of logC*^'^'^^^ can be used in the global fit to improve the accuracy of the 
determination of R123, which in this scenario turns to be 

^mi:eTiar = 1-004 ±0.007. (4.19) 

Here (/A'//7r)//+(0) has been fixed from lattice. This ratio is the key quantity to be 
improved in order to reduce present uncertainty on Ri23- 

The measurement of R123 above can be used to set bounds on the charged Higgs mass 



and tan /3. Figure IC shows the excluded region at 95% CL in the M/^- tan /3 plane (setting 
eo = 0.01). The measurement of BR(i? tu) ||7^ can be also used to set a similar bound 
in the M//-tan /3 plane. While B ^ tv can exclude quite an extensive region of this plane, 
there is an uncovered region in the exclusion corresponding to a destructive interference 
between the charged-Higgs and the SM amplitude. This region is fully covered by the 
K ^ result. 



In the case of right-handed currents [12|, R123 can be obtained from a global fit to 



the values of eqs. ( [4.1| ) and ( [4. 2D . Here logC^^P is free of new physics effects and can be 
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also used to constrain ifK/fn)/f+{^) together with lattice results (namely the values in 



tab. ^0|). The result is 

^;23lXurr. = 1-004 ± 0.006 . (4.20) 

4.4 Tests of Lepton Flavor Universality 

4.4.1 Lepton universality in K^^ decays 

The test of Lepton Flavor Universality (LFU) between K^^ and JT^a modes constraints a 
possible anomalous lepton-flavor dependence in the leading weak vector current. It can 
therefore be compared to similar tests in r decays, but is different from the LFU tests in 
the helicity-suppressed modes 7r/2 and K12. 

The results on the parameter r^e = ^^^^/k sZ-^^s/k 3 

r^e = 1-004 ± 0.004 , (4.21) 

in excellent agreement with lepton universality. Furthermore, with a precision of 0.5% the 
test in decays has now reached the sensitivity of r decays. 

4.4.2 Lepton universality tests in K£2 decays 

The ratio Rk = r(i^^2)/r(-f^e2) can be precisely calculated within the Standard Model. 
Neglecting radiative corrections, it is given by 

< = 4 = 2.569 X 10-^ (4.22) 



and reflects the strong helicity suppression of the electron channel. Radiative corrections 
have been computed with effective theories B, yielding the final SM prediction 



Rf^ = i?5J)(l+5i?5^d.corr. 

2.569 X 10"^ X (0.9622 ± 0.0004) = (2.477 ± 0.001) x 10"^ . (4.23) 



Because of the helicity suppression within then SM, the Ke2 amplitude is a prominent 
candidate for possible sizable contributions from physics beyond the SM. Moreover, when 
normalizing to the ir^2 rate, we obtain an extremely precise prediction of the Ke2 width 
within the SM. In order to be visible in the Ke2/K^2 ratio, the new physics must violate 
lepton flavor universality. 

Recently it has been pointed out that in a supersymmetric framework sizable viola- 
tions of lepton universality can be expected in K12 decays [^]. At the tree level, lepton 
flavor violating terms are forbidden in the MSSM. However, these appear at the one-loop 
level, where an effective H^luj- Yukawa interaction is generated. Following the notation of 



Ref. (see also Section 22), the non-SM contribution to Rk can be written as 

(4.24) 



Rk ~ Rr 



The lepton flavor violating coupling A13, being generated at the loop level, could reach 
values of 0(10"^). For moderately large tan/? values, this contribution may therefore 
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Figure 11: Exclusion limits at 95% CL on tan/3 and the charged Higgs mass Mh± from 
\Vus\Ki2/\Vus\Kn for different values of A13. 



enhance Rk by up to a few percent. Since the additional term in Eq. [4.24| goes with the 
forth power of the meson mass, no similar effect is expected in 1^12 decays. 

The world average result for Rk presented in Section |3| gives strong constraints for 
tan (3 and Mu± , as shown in Fig. For values of A13 ~ 5 x 10"'^ and tan /3 > 50 the 
charged Higgs masses is pushed above 1000 GeV/c^ at 95% CL. 
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A. BRS fit procedure 



The fits to Kl and data are performed with FORTRAN programs. MIGRAD is used for 
the minimization; errors are obtained with MINOS. 

Suppose we have A'^ measurements of M quantities, e.g., BRs, ratios of BRs, Hfetimes, 
or partial widths, where N > M as some quantities arc measured by more than one 
experiment. Denote the A'' measurements Xj, and the expected value for each as calculated 
from the free parameters of the fit Xj. We also refer to the expected values for quantities 
measured by more than one experiment by the index m, i.e., Xm with m = 1, M. 

The errors on the inpiit parameters arc denoted o"j. All errors on the input parameters 
are assumed to be Gaussian. For uncorrelated measurements with statistical and systematic 
errors quoted separately, we add the errors in quadrature. In many cases, the results for 
different quantities measured by the same experiment have correlated errors. The errors 
are then described by the covariance matrix Vij, with Va = af and Vij = pijaiaj. The 
expression to be minimized is then 

N N 

X' = J2 E(^^ - ^^)(^^- - ■ (A.1) 

i=l j=l 

In practice, V is block diagonal and only the relevant sub-matrices are inverted. 

The penalty method is used to implement the constraint on the sum of the BRs. In this 
method, a term G (1 — ^ BR)^ is added to the to be minimized. As G is increased, the 
constraint is enforced with greater and greater precision and the result of the fit saturates 
(until at some very large value of G, problems related to the precision of the calculation 
set in). G is determined by trial and error; its value is 2 X 10^ for the Kl fit and 1 x 10^ 
for the fit. The Kl fit is somewhat more sensitive to the value of G, because the Kl 
BRs entering the fit span three orders of magnitude. As a result, precision problems have 
a greater effect on the constraint balance. 

Once the fit has been performed, scale factors are calculated and used as described in 
the general introduction to the PDG compilation. As above, our N data points consist 
of m = 1, M distinct measured quantities, each of which is measured by rim experiments, 
indexed by k^,- {N = n^). Here it is useful to adopt the notation Xmkm^'^mkm ^^r the 
individual measurements, and the notation Xm for the expected value for the m**^ quantity. 
After the fit is performed once, the error a^n on Xm is evaluated from the output covariance 
matrix for the free fit parameters. Then, the scale factor for the measured quantity m is 
calculated as 

-i / _ \2 

q2 _ ^ \ " {^mkm ~ ^m) / « r)\ 

Next, the errors (Tmkm. scaled by the greater value of Sm and unity. For subsets of 
correlated measurements (all from the same experiment), the index km can be omitted 

to write V„im' = Prnrn'^mf^rn'] thc scalc factors are applied to o"m and o"„j' and V and its 
inverse are recalculated. Finally, thc fit is performed a second time. For each of the fit 
parameters, we report the central value from the first fit, and the error (and correlations) 
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from the second fit. The scale factors for the errors on the fit parameters are defined as the 
ratios of the errors from the second fit to those from the first. The value of reported is 
from the first fit. 

For the purposes of comparison, pull values are calculated for each measurement simply 
as {xi - Xi)/ai. 

For the BR/lifetime fits, the errors are in general symmetric to within rounding error; 
in any case we report the greater of the positive and negative minos errors. 



B. Fit for BRs and lifetime 

The 8 free parameters in the Kl fit are BR{Ke3), BR(i^^3), BR(37rO), BR(7r+7r-7r°), 
BR(7r"'"7r~), BR(7r''7r''), BR(77), and tk^- The fit makes use of the 18 measurements 
in Table |ll|. With one constraint, the fit has 11 degrees of freedom. 

The differences between our fit and the 2006 PDG fit are as follows: 

• In our fit, the intermediate KTeV and KLOE values (i.e., before applying constraints) 
are the inputs, and the complete error matrix is used to handle the correlations 
between the measurements from each experiment. In the 2006 PDG fit, the final 
KTeV and KLOE BR results were used and one measurement involving BR(37r'') 
was removed in each case. 



Our fit makes use of the preliminary BR(37r*^) 
measurements from NA48. 



and new BR(7r+7r-)/BR(if, 



e3j 



Our fit parameter BR(7r"'~7r~) is understood to be inclusive of the DE component. 
This helps to satisfy the constraint. The input data are treated consistently in this 
respect. 

We do not make use of the measurement of BR(77) /BR(7r'^7r'^) from NA31 (Burkhardt 
'87), since both we and the PDG have excluded the other measurements from NA31. 



Notes on data in Table 11 



1. Direct measurement of r/^^ from Sir^ events; independent of other KLOE measure- 



ments 1 31]. 



2. We make use of the KLOE results for the main Kl BRs (#3, #4, #6, and #9) 
obtained before applying constraints on the sum of the BRs |3^]. The BR values in 
Table ^ thus depend on tk^ as follows: 



BR 



BR' 







1 + 0.0128 ns-i r' 



Kl 



where 



51.54 ns. The errors listed for these values in Ref. [^] include an explicit 
contribution from the uncertainty on the reference value of tk^ ■ This contribution has 



been unfolded from the errors in Table 11. In addition, these four BR measurements 
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Point Parameter 



Value 



Source 



Note 



1 ^Kl 

2 TK^ 

3 BR(i^e3) 

4 BR(K^3) 

5 BR{K^3)/BR{K,s) 

6 BR(37rO) 

7 BR(37rO)/rK^ 

8 BR(37rO)/BR(is:e3) 

9 BR(7r+7r-7rO) 

10 BR(7r+7r-7rO)/BR(i^e3) 

11 BR(7r+7r-)/BR(ii:e3) 

12 BR(7r+7r-)/BR(i^e3) 

13 BR(7r+^-)/BR(i^^3) 

14 BR(i^e3)/BR(2 tracks) 

15 BR(7r07rO)/BR(37rO) 

16 BR(7r07rO)/BR(7r+7r-) 

17 BR(77)/BR(37rO) 

18 BR(77)/BR(37rO) 



50.92(30) ns 
51.54(44) ns 
0.4049(21) 
0.2726(16) 
0.6640(26) 
0.2018(24) 
3.795(58) MHz 
0.4782(55) 
0.1276(15) 
0.3078(18) 
0.004856(29) 
0.004826(27) 
0.007275(68) 

0.4978(35) 
0.004446(25) 
0.4391(13) 
0.00279(3) 
0.00281(2) 



KLOE '05 
Vosburgh '72 
KLOE '06 
KLOE '06 
KTeV '04 
KLOE '06 
NA48 '04 
KTeV '04 
KLOE '06 
KTeV '04 
KTeV '04 
NA48 '06 
KLOE '06 
NA48 '04 
KTeV '04 
PDG '06 
KLOE '03 
NA48 '03 



1 m 
m 



2 Igi 



2 g 



3 12| 



2 g 



4 m 

3 H 
2 ISOtl 



3 l2^ 
3,5 [|7 
5 PH 



5 PI 



6 |28|1 



3 |27|] 



7 m 



m 



Table 11: Input data used for the fit to BRs and lifetime. 

are correlated by common systematics as described in KLOE Note 204, although the 
full correlation matrix is not given therein. The correlation matrix is as follows: 

3 4 6 9 

3 1.000 0.091 0.069 0.494 

4 1.000 -0.025 0.267 
6 1.000 0.074 
9 1.000 

3. The correlation matrix for the KTeV relative BR measurements (#5, 7^8, #10, #11, 



and #15) is as follows [27|: 



5 8 10 11 15 

5 1.00 0.14 0.21 0.24 0.09 

8 1.00 -0.06 -0.07 0.30 

10 1.00 0.49 0.04 

11 1.00 0.07 
15 1.00 



4. This is based on the preliminary NA48 measurement BR(37r'') = 0.1966(34), as re- 
ported in Ref. [29|. R. Wanke has confirmed that the 2004 PDG value for tkj^ was 
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used to obtain this result. The NA48 value for this BR scales directly with the life- 
time value used. R. Wanke has supplied the value in the table for the partial width, 
with the contribution to the error on the BR from the uncertainty on the lifetime 
unfolded. 

5. The fit value of BR(7r+7r^) includes the DE component. The KLOE measurement of 
BR(7r+7r-)/BR(K^3) |3| (#13) is inclusive of DE. The KTeV and NA48 measure- 
ments of BR(7r+7r^)/BR(i^e3) (#11 and #12, respectively) are treated as follows. 



We use the average values of DE/(DE + IB) from Refs. [SQ] and Q together 
with BR(7r+7r-7iB,£^7 > 20 MeV)/BR(7r+7r-) = 7.00 x 10"^ ||, to calculate 
that DE accounts for 1.52(7)% of the inclusive Kl tt^tt^ width. The error 
on this correction is negligible for the purposes of the fit. 

The KTeV measurement of BR(7r+7r~)/BR(/i:e3) (#11) excludes DE (in the 



sense that Ref. |27| says that DE is not in the generator for the acceptance 
calculation). We therefore subtract 1.52% from the fit value of the ratio when 
calculating the contribution to from this KTeV measurement. 

• The contribution from DE to the NA48 measurement BR(7r^7r^)/BR(A'e3) = 
4.835(22) (20) x 10^^ is estimated to be 0.19(1)%, which we subtract to obtain 
value #12. We then handle the data point in the same way as we do for KTeV. 

6. For our fit, BR(2 tracks) has to be calculated from the free fit parameters. Like the 
PDG, we use 

BR(2 tracks) = BR(Ke3) + BR(K^3) + 0.03508 BR(37r°) 
+ BR(7r+7r"7r°) + BR(7r+7r"). 

7. From the ETAFIT analysis |4|]. 
B.l Results 



The results of the fit are summarized in Table 12. The output correlation matrix is given in 



Table 13. The pull values for the input measurements are listed in Table 0. With respect 
to the 2006 PDG fit, our fit has a somewhat lower probability. 

When our fit is run without inclusion of points #7 and #12, without DE corrections 
for the vr^vr^ channel, and with the measurement of BR(77)/BR(7r'^7r'^) from NA31, we 
reproduce the 2006 PDG fit result. In this configuration, the only difference between our 
fit and the 2006 PDG fit is the treatment of the BR and lifetime data from KLOE and 
KTeV. We obtain the same values for all 8 fit parameters, with x^/ndf = 14.9/10. Our 
scale factors in this case are more uniform than those obtained in the 2006 PDG fit; in 
particular, for BR(i^e3), BR(37rO), and BR(7r+7r~7r°) we have S = 1.2, 1.1, and 1.4, to be 
compared with the second column of Table p^ . 

Excluding the measurement of BR(77)/BR(7r^7r'^) from NA31 has a negligible effect on 
the fit results, while the number of degrees of freedom is reduced by one, giving x^/ndf = 
14.9/9 (9.4%). Turning on the DE correction degrades the fit quality from x^/ndf = 14.9/9 
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This fit 2006 PDG 

18 measurements 17 measurements 



xVndf = 19.7/11 (4.9%) xVndf = 14.8/10 (14.0%) 



Parameter 


Result 


S 


Result 


S 


BR(Ke3) 


0.4058(9) 


1.3 


0.4053(15) 


2.1 


BR(K^3) 


0.2706(8) 


1.3 


0.2702(7) 


1.0 


BR(37r°) 


0.1943(10) 


1.3 


0.1956(14) 


1.9 


BR(7r+7r-7rO) 


0.1259(8) 


1.5 


0.1256(5) 


1.0 


BR(7r+7r-) 


1.986(7) X 10-3 


1.2 


1.976(8) X 10-3 


1.0 


BR(7r07rO) 


8.60(5) X 10-^ 


1.5 


8.69(4) X 10-"^ 


1.1 


BR(77) 


5.45(4) X 10"^ 


1.1 


5.48(5) X 10-"^ 


1.2 


BR(rA-J 


51.15(20) ns 


1.1 


51.14(21) ns 


1.0 



Table 12: Results of fit to Kl BRs and lifetime, with comparison to 2006 PDG fit. 



-0.286 -0.422 


-0.288 


+0.112 


-0.282 


-0.270 


-0.005 


+1.000 -0.378 


-0.217 


-0.046 


-0.216 


-0.241 


+0.183 


+1.000 


-0.354 


-0.029 


+0.609 


+0.637 


-0.036 




+1.000 


-0.035 


-0.205 


-0.226 


-0.127 






+1.000 


+0.205 


-0.020 


-0.033 








+1.000 


+0.387 


-0.029 










+1.000 


-0.027 












+1.000 



Table 13: Correlation matrix for output parameters of Kl fit. 



to 19.6/9 (2.02%). When points #7 and #12 are added, the fit quality is slightly improved 
and the result in Table |l^ is obtained. 

However, the fit quality improves dramatically when the PDG ETAFIT result for 
BR(7r'^7r'^)/BR(7r+7r^) (#16) is removed. This is true independently of whether or not the 
DE correction and/or the additional results are included. For example, our same fit with 
the PDG ETAFIT point removed gives x^/ndf = 14.8/10 (13.9%), with changes in the fit 
values for the BRs at the la level. In all other configurations (DE correction on/off; points 
#7, #12, NA31 BR(77)/BR(7r°7rO) included/excluded), the fit gives similar resuhs. 

Using the values of BR(7r+7r-) and BR(7r°7r°) fr om our fits including and excluding 
the PDG ETAFIT point, we have evaluated Re e'/e from 



Re e'/e = - 
' 6 



BR(^°7r°) 
BR(7r+7r-) 



1-^5 



with Rs = BR(Ks 7r+7r-)/BR(Ks ^ vr°7r°) = 2.2549(54) [||] as described in Sec. U 
We obtain 

Re e'/e = (-25 ± 23) x IQ-^ (without ETAFIT); 

Re e'/e = (14 ± 11) x 10-^ (with ETAFIT); 
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Point 


Parameter 


Source 


Pull 


6 


BR(3^0) 


KLOE '06 


+2.74 


2 




Vosburgh '72 


+0.88 


15 


BR(7r0^0)/BR(37r0) 


KTeV '04 


+0.81 


9 


BR(7r+7r-7r°) 


KLOE '06 


+0.71 


4 


BR(i^^3) 


KLOE '06 


+0.41 


18 


BR(77)/BR(3^0) 


NA48 '03 


+0.31 


12 


BR(7r+^-)/BR(Ke3) 


NA48 '06 


+0.22 


7 


BR(3^0)/rKL 


NA48 '04 


-0.07 


8 


BR(3^0)/BR(i^e3) 


KTeV '04 


-0.13 


17 


BR(77)/BR(37r°) 


KLOE '03 


-0.46 


16 


BR(7r07r°)/BR(7r+7r-) 


PDG '06 


-0.57 


14 


BR(i<"e3)/BR(2 tracks) 


NA48 '04 


-0.71 


1 




KLOE '05 


-0.78 


13 


BR(7r+^-)/BR(K^3) 


KLOE '06 


-0.94 


5 


BR(K^3)/BR(ire3) 


KTeV '04 


-1.11 


11 


BR(7r+7r-)/BR(Ke3) 


KTeV '04 


-1.32 


3 


BR(i^e3) 


KLOE '06 


-1.37 


10 


BR(7r+7r-7rO)/BR(ire3) 


KTeV '04 


-1.39 



Table 14: Pull values for input data used in fit to BRs and lifetime. 



to be compared to the current PDG average, (16.7+2.3) xlO""^. The ETAFIT point IS very 
precise; when it is included, the fit results for BR(7r'^7r'^)/BR(7r'^7r~) are highly constrained. 
This pulls down the value of BR(7r'^7r'^), and, also of BR(37r'^), via the KTeV measurement 
of BR(7r'^7r'^)/BR(37r'^). As seen from Table 14, the measurement with the largest positive 
pull on the fit is the KLOE measurement of BR(37r'^), which PDG has chosen to eliminate 
from the 2006 fit as part of their treatment of the correlated KLOE measurements. 

We emphasize that the values of BR(iire3) and BR(if^3) are not affected very much 
by these developments. The scale factor, and hence the reported error, on BR(i^e3) is 
significantly smaller in our fit, which spreads the pulls somewhat more evenly over the 
different measurements than does the PDG fit. 



C. Fit for K± BRs and lifetime 

The 7 free parameters in the fit are BR(i^^2), BR(7r7rO), BR(7r7r7r), BR(ii'e3), BR(K^3), 
BR(7r7r°7r'^), and Tx±- The fit makes use of the 26 measurements in Table With one 
constraint, the fit has 24 degrees of freedom. 

The principal difference between the fit performed here and the 2006 PDG fit is that 
our fit includes the following recent measurements: 

• Preliminary tx± from KLOE (7^6, #7); 

• Preliminary BR{Kez) and BR(K/,3) from KLOE (#14, #20); 
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Point Parameter 



Value 



Source 



Note 



1 
2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 



Tk± 
Tk± 
Tk± 
Tk± 
Tk± 
Tk± 

BR(7r7r' 
BR(7r7r' 
BR(7r7r' 
BR(7r7r' 
r(7r7r7r) 

BR(i^e3 
BR(i^e3 

BR(i^e3 

BR(i^;.3 
BR(i^;.3 
BR(i^^3 
BR(i^^3 
BR(i^^3 
BR(7r7r°7r°) 
BR(7r7r07r0)/BR(7r7r7r 



/BR(K^2) 
/BR(K^2) 
/BR(K^2) 



/BR(7r^o+K^3+7r27rO) 
/BR(K^2 + ^vrO) 
/BR(K^2 + ^vrO) 
/BR(7r7rO) 
/BR(7r7rO) 

/BR(7r7r°) 

/BR(Ke3) 
/BR(Ke3) 
/BR(Ke3) 



12.451(30) ns 
12.368(41) ns 
12.380(16) ns 
12.272(36) ns 
12.443(38) ns 
12.367(78) ns 
12.391(55) ns 
0.6366(17) 
0.2066(11) 
0.3329(48) 
0.3355(57) 
0.3277(65) 
4.513(24) MHz 
0.04965(53) 
0.1962(36) 
0.0616(22) 
0.0589(21) 
0.2449(16) 
0.2470(10) 
0.03233(39) 
0.1636(7) 
0.671(11) 
0.670(14) 
0.667(17) 
0.01763(26) 
0.303(9) 



Koptev '95 [33 1 



Koptev '95 [33 



Ott '71 |4[ 
Lobkowicz '69 (H 
Fitch '65 M 



KLOE '06 g| 



KLOE '07 H 



KLOE '06 |g 
KLOE '07 gl| 
Usher '92 [|]] 
Weissenberg '76|8^ 
Auerbach '67 ||] 
Ford '70 



KLOE '07 |39|| 
Sher '03 
Eschstruth '68 
Cester '66 [^ 
ISTRA '07 



NA48 '07 |3[ 
KLOE '07 pi 



NA48 '07 Q 
Horie '01 [||] 
Heintze '77 [p 



Botterill '68 [g6[ 



KLOE '04 H 



Bisi '65 |8[ 



2,4 



2 

5 

2,4 
5 



Table 15: Input data used for the fit to K BRs and lifetime. 

• Preliminary BR(7r7r°) from KLOE (#9); 

• Preliminary BR(i^e3)/BR(7r7r°) from ISTRA+ (#18); 

• BR(Ke3)/BR(7r7r°) (#19) and BR(K^3)/BR(7r^0) from NA48/2 (#21). 

These new measurements have a profound impact on the results of the fit. Other differences 
are as follows. 

• In the 2006 PDG fit, BR(7r''7r''ei^) is a free parameter (but curiously, BR(7r7rei^), for 
which there is a published measurement from E865 with much higher accuracy [p9|[ , 
is not). The PDG fit therefore uses three measurements involving BR(7r'^7r'^ei^) and 
BR(7r''7r'^ez^)/BR(iire3) that are not used in our fit. 
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• We don't use the six BR measurements from Chiang '72. Our reading of Chiang 
'72 suggests that no attempt was made to implement radiative corrections for the 
branching ratio analysis. In addition, the six BR measurements from Chiang '72 are 
constrained to sum to unity. The correlation matrix is not available. PDG omits 
BR(7r7r7r). 

It would be highly desirable to discard many other old measurements in the fit as 
2006 PDG has done for the Kl fit. Unfortunately, are no recent measurements involving 
BR(7r7r7r). As a result, the fit is unstable if only recent measurements are used. 



Notes on data in Table 15 



1. The only difference between the Koptev measurements is the material used for the 
kaon stopper (#1-U, #2-Cu). 

2. Preliminary measurement. 

3. The dependence of these BRs on the lifetime is accounted for in the fit: 

BR = BR° [1 + 0.0405(rj^± - r]^±)] 

where BR*' is evaluated with t^± = 12.360 ns. The uncertainty from the value of 
Tj^± may not have been properly unfolded. In addition, these two measurements 
are have a correlation coefficient of 0.627, mainly from the use of common efficiency 
corrections. 

4. The recent NA48 pubhcation |3^ gives values for BR(iCe3)/BR(7r7rO), 
BR(K^3)/BR(7r7rO). The value of BR(i^e3)/BR(7r7r°) has been updated at KAON07. 

C.l Results 



The results of the fit are summarized in Table 16. The output correlation matrix is given 



in Table 17. The pull values for the input measurements are listed in Table The poor 
fit quality derives from the following sources. 

• The fit quality is significantly degraded by the scatter in the five older measurements 
of T^±; when these are replaced with their PDG average with scaled error, t^± = 
12.385(25) ns, the fit gives x^/ndf = 24.3/16 (8.4%), with no significant changes in 
the results. Note that after this treatment the fit quality is about the same as it is 
for the 2006 PDG fit (which, however, includes all of the older r^± measurements 
without taking the average). 

• There is some confiict among the newer measurements involving BR(i^e3)! ^ seen 
from the pulls for the NA48 '07 (#19), Sher '03 (#15), ISTRA '07 (#20), and KLOE 
'07 (#14) measurements: +1.04, -0.26, -0.74, and -2.13, respectively. (Table |l8|) . 



The evolution of the average values of the BRs for decays and for the important 
normalization channels as a result of the introduction of the preliminary measurements is 
evident in Fig. |2[ The figure dramatically illustrates why experiments that measure ratios 
such as BR(i^e3)/BR(7r7r*') should always quote the ratio with usable errors, in addition to 
the normalized, final value. 
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This fit 
26 measurements 



2006 PDG 
26 measurements 





xVndf = 42/20 (0.31%) 


xVndf = 30/19 (5.2%) 


Parameter 


Result 


c 
D 


Result 


c 
D 




63.57(11)% 


1.1 


63.44(14)% 


1.2 


BR(wO) 


20.64(8)% 


1.1 


20.92(12)% 


1.1 


BR(7r-7r7r) 


5.595(31)% 


1.0 


5.590(31)% 


1.1 


BR(i^e3) 


5.078(25)% 


1.2 


4.98(7)% 


1.3 




3.365(27)% 


1.7 


3.32(6)% 


1.2 


BR(7r7r°7rO) 


1.750(24)% 


1.1 


1.757(24)% 


1.1 


BR(7r°7rOez^) 


Not in fit 




2.2(4) X 10"^ 


1.0 


BR(t^±) 


12.384(19) ns 


1.7 


12.385(24) ns 


2.1 



Table 16: Results of fit to BRs and lifetime, with comparison to 2006 PDG ht. 



1.000 



-0.874 
1.000 



-0.170 
-0.121 
1.000 



-0.725 


-0.548 


-0.258 


-0.045 


0.610 


0.333 


0.031 


-0.032 


-0.100 


-0.074 


0.055 


0.273 


1.000 


0.442 


0.009 


-0.030 




1.000 


-0.010 


-0.020 






1.000 


0.010 








1.000 



Table 17: Correlation matrix for output parameters of fit. 



D. Averages of form- factor slopes 
D.l Procedure 

We work principally with quadratic form- factor slope parametrization. To average the 
form-factor slopes, a fit with correlations is performed. Scale factors for the errors are 
calculated as described in section For the fit to the form-factor slopes, since there are 
no measurements of combinations of the fit parameters, the scale factors can be obtained 
directly from Eq. ( [A.2| ). Because of the high degree of correlation in the measurements 
of A' and A", a large scale factor may result in a small change in from the fits. We 
therefore report scaled errors only when the value of x^/ndf is unsatisfactory. 



D.2 Input data 



The data used in the fit are summarized in Table |19|. The following notes apply to the 
table entries. 



1. In our combined fits to i^es and iC^s data, we use the averages quoted by KLOE and 
KTeV rather than using their K^^ and K^^ measurements separately. In any event, 
our averages of the i^es and K^^ results from each experiment have good values of 
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Point Parameter 



Source 



Pull 



1 




Koptev '95 


+2.25 


11 


BR(7r7rO)/BR(i^^2) 


Weissenberg '76 


+1.89 


10 


BR(7r7rO)/BR(i^^2) 


Usher '92 


+1.70 


5 




Fitch '65 


+1.56 


21 


BR(K^3)/BR(7r7rO) 


NA48 '07 


+1.04 


19 


BR(Ke3)/BR(7r7rO) 


NA48 '07 


+1.03 


22 


BR(K^3)/BR(Ke3) 


Horie '01 


+0.76 


16 


BR(Ke3)/BR(i^^2 + vtttO) 


Eschstruth '68 


+0.59 


8 


BR(K;.2) 


KLOE '06 


+0.52 


23 


BR(K^3)/BR(Ke3) 


Heintze '77 


+0.52 


25 


BR(7r7r°7r°) 


KLOE '04 


+0.52 


12 


BR(7r^0)/BR(/C^2) 


Auerbach '67 


+0.46 


24 


BR(K^3)/BR(Ke3) 


Botterill '68 


+0.26 


7 




KLOE '07 


+0.14 


13 


r(7r7r7r) 


Ford '70 


-0.22 


6 




KLOE '06 


-0.21 


3 




Ott '71 


-0.22 


15 


BR(Ke3)/BR(i^^3 +vr7r0 +7r27r0) 


Sher '03 


-0.26 


2 




Koptev '95 


-0.38 


17 


BR(Ke3)/BR(i^^2 + ^vrO) 


Cester '66 


-0.67 


20 


BR(Ke3)/BR(7r7rO) 


ISTRA '07 


-0.74 


26 


BR(7r7r°7r°)/BR(7r7r7r) 


Bisi '65 


-1.07 


14 


BR(Ke3) 


KLOE '07 


-2.13 


4 




Lobkowicz '71 


-3.10 


20 


BR(K^3) 


KLOE '07 


-3.41 



Table 18: Pull values for input data used in fit to BRs and lifetime. 



X^/ndf and confirm the results quoted by the experiments, including the correlation 
coefficients. 

2. The exact value of p{\'j^, A'^) is not available for the NA48 K^z measurement. NA48 
and PDG together estimated p = —0.88; this value appears in the 2006 PDG listings 
Q. For use with Eq. (p|), we put A'| = 2)^'}^^^^\ 

3. An official value of p(AY,A") is not available for the ISTRA+ K^z measurement; 
the value in the Table was obtained directly from the collaboration. For use with 
Eq. (pD, we put A^ = C aJ™^^ and A'| = 2C2A?™^\ with C = {m^+/m^of = 
1.069223. 

4. Systematic errors for the ISTRA+ quadratic fit results for i^^3 are not given in 
Ref . |5^] ; the errors in the table are statistical only. Nor are the correlation coefficients 
available; these have been obtained directly from the collaboration. For use with 
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Eq. (2^), A'_)_ and A'_{_ are converted as above; we also put Xq = C Xq^"^^^K Finally, 
we note that no information concerning the treatment of radiative corrections is 
given in Ref. ||5^ . Failure to account for radiative effects could result in a noticeable 
systematic shift in the slope results. 



For the KLOE and KTeV form-factor slope measurements, the correlation coefficients 
apply to the total errors (statistical and systematic). For the ISTRA-I- and NA48 K^^ 
slopes, the correlation coefficients appear to apply to the statistical errors. In our fits, we 
assume that the correlation coefficients apply to the total errors on the form-factor slopes 
(statistical and systematic). This approximation is motivated as follows. In general, the 
systematic errors are estimated by varying analysis parameters and refitting. In that case, 
the statistical correlations naturally present will also affect the excursions due to systematic 
variations, see Appendix ^ 



D.3 Fit results for K13 slopes excluding NA48 Kfj,s data 



Kl and K 



Kl only 



The result of our fit to all data is presented in Table 0. As discussed in Sec. 3.5.2 , the 
NA48 form-factor slope measurements are in contrast with the results from the other 
experiments. As an exercise, we fit all results in Table 19 except the NA48 measurement of 
the slopes |^^. The results are shown in Table 2C. The first column of the table gives 
the results of the fit to all other measurements from KLOE; the second gives the results 
of the fit to the Kl measurements from KLOE, KTeV, and the Kl es measurement from 
NA48. 

The evaluations of the phase-space 
integrals for all four modes are listed 
in each case. Correlations are fully ac- 
counted for, both in the fits and in the 
evaluation of the integrals. The val- 
ues of x^/ndf do not raise any signifi- 
cant concerns about the compatibility 
of the input data. The fit to all data 
gives xVndf = 12.6/10 (25.0%). 

The evaluations of the phase-space 
integrals for all four modes are listed 
in each case. Correlations are fully ac- 
counted for, both in the fits and in the 
evaluation of the integrals. The val- 
ues of x2/ndf do not raise any signifi- 
cant concerns about the compatibility 
of the input data. The fit to all data 
gives x2/ndf = 12/9 (22.3%). 



Measurements 
xVndf 
A'^ X 10^ 
A'l X 10^ 
Ao X 10^ 
P(A'+,A'|) 
p(AV,Ao) 
p(A'|,Ao) 

ml,) 



13 

13/9 (24.9%) 
25.0 ±0.8 
1.6 ±0.4 
16.0 ±0.8 
-0.94 
±0.26 
-0.37 
0.15459(20) 
0.15894(21) 
0.10268(20) 
0.10559(20) 
±0.59 



9/5 (12.3%) 
24.5 ± 1.1 
1.8 ±0.4 
14.8 ± 1.1 
-0.95 
±0.28 
-0.38 
0.15446(27) 
0.15881(28) 
0.10236(28) 
0.10532(29) 
±0.62 



Table 20: Averages of quadratic fit results for Ke3 and 
Kfj,3 slopes, excluding new K^^ data from NA48. 
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E. Error estimates 



It is quite easy to estimate the ideal error in the measurements of a set of parameters 
p={pi, P2t ■ ■ Pn) from fitting some distribution function to experimentally determined 
spectra. Let -F(p, x) be a probability density function, PDF, where p is some parameter 
vector, which we want to determine and x is a running variable, like t. The inverse of the 
covariance matrix for the maximum likelihood estimate of the parameters is given by [100| : 

'■^ dpidpj 
from which, for N events, it trivially follows: 



F dpi dp 



with dv the appropriate volume element. We use in the following the above relation to 
estimate the errors on the FF parameters for one and two parameters expression of the FFs 
f+{t) and fo{t). The errors in any realistic experiment will be larger than our estimates, 
typically two to three times. The above estimates are useful for the understanding of the 
problems in the determination of the parameters in question. The elements of G depend on 
the values of the parameters. In the case of the form factors, the errors on the A parameters 
change insignificantly for 10% changes of the parameters. In other words the errors do not 
depend on the data being fitted and the correlations apply also to the systematic part of 
the errors. 

E.l Ke3 decays 

For a quadratic FF, f{t) = 1 + \'_^{t/m?') + {\'^/2){t/m'^)'^, the inverse of the covariance 
matrix G^^, the covariance matrix G4. and the correlation matrix are: 

/ 5.937 13.867 \ 1 / 1.258^ -0.606 \ / 1 - 
\ 13.867 36.2405/ -0.606 0.509^ J ' ^ 

The square root of the diagonal elements of G+ gives the errors, which for one million 
events are (5A'_,_=0. 00126, =0.00051. The correlation is very close to —1, meaning that, 
because of statistical fluctuation of the bin counts, a fit will trade A'_|_ for A" and that the 
errors are enlarged. A fit for a linear FF, f{t) = 1 + X'j^{t/m?) in fact gives A^=0.029 
instead of 0.025 and an error smaller by ~3: 

5AV = VG+(1,1) = 0.0004. 

A simple rule of thumb is that ignoring a term, increases A^ by ~3.5x A'^. For Kf.'i decays 
the presence of a term in the FF is firmly established. It is however not fully justified 
to fit for two parameters connected by the simple relation A^=2xAY^. The authors of ref. 
explicitly give an error for their estimate of the coefficient of the terms. 
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E.2 K^-i decays 

The scalar FF only contributes to i^^s decays. Dealing with these decays is much harder 
because: a) - the branching ratio is smaller, resulting in reduced statistics, b) - the E^^ or t 
range in the decay is smaller, c) - it is in general harder to obtain an undistorted spectrum 
and d) - more parameters are necessary. This is quite well evidenced by the wide range of 
answers obtained by different experiments p8|. 



54, p3]. Assuming that both scalar and 



vector FF are given by quadratic polynomials as in Eq. ( |2.5D , ordering the parameters as 
A'o, Aq, and A" , the matrices Gq^_(_ and are: 



N 



/1.64 
5.44 
1.01 

V3.90 



5.44 
18.2 
3.01 
12.3 



1.01 
3.01 
1.47 
4.24 



3.90 \ 
12.3 
4.24 
13.8/ 



1 



/ 63.9^ 
1200 
-923 
V 197 



-1200 
18.82 
272 

-59 



-923 
272 
14.8^ 

-48 



197 \ 
-59 
-49 
3.4V 



and the correlations, ignoring the diagonal terms, are: 

/ -0.9996 -0.974 0.91 



0.978 



V 



-0.919 
-0.976, 



(E.l) 



All correlations are very close to —1. In particular the correlations between Aq and Aq is 
—99.96%, reflecting in vary large JAq and 5\'q errors. We might ask what the error on 
Aq and Aq might be if we had perfect knowledge of A^ and A" . The inverse covariance 
matrix is give by the elements (1,1), (1,2), (2,1) and (2,2) of the Gq]^^ matrix above. The 
covariance matrix therefore is : 



Go(A^, K for AV, A'l known) = - 




-20 
2.42 



For one million events we have (5Aq=0.0024, about 4x the expected value of Ag. In other 
words Aq is likely to be never measurable. It is however a mistake to assume a scalar 
FF linear in t, because the coefficient of t will absorb the coefficient of a term, again 
multiplied by ~3.5. Thus a real value Ag=0.014 is shifted by the fit to 0.017, having used 



the parametrization of Ref . [ 12 1 . 



E.3 From the linear to the dispersive parametrization 

The results on the FFs obtained with the linear parametrization can be used to determine 
the parameter of the dispersive parametrization. As shown in the previous section the 
correlation between Ag and Aq is close to -1 and any fit to Ag, Aq form K^-^ decays will give 
values satisfying the relation: 

a;,' = tan^AJ) + 5; (E.2) 
tan (p is independent on the number of events of the experiment 

2pi,20'2cri 



tan 2(f) 



da 



with 1 = A^, 2 = X'^ 



(E.3) 
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0.002 



0.001 




0.02 



Figure 12: Linear parametrization extrapolation along correlation line and relation from dispersive 
parametrization. 



and amounts to -0.3. B can be determined from the experimental results for Aq obtained 
using the linear parametrization ( Aq=0). 

Therefore we can translate the results for Aq in any new parametrization with only one 
parameter with almost negligible 3th order term. In particular the dispersive parametriza- 
tion gives: 

A[)' = A^V (4.16 ± 0.50) X 10"^ (E.4) 
the procedure is shown in figure 12 for Ao=(15 it 1) x 10^'^. 



49 



O 



Experiment 


AV X 10^ 


A'l X 10^ 


Ao X 10^ 


P(AV,A'|) 


p(AV,Ao) 


/'(A'|,Ao) 


Analysis 


Note 


KLOE Kl e3 [49] 


25.5 ± 1.8 


1.4 ±0.8 




-0.95 






t from Ks 7r"^7r~ 




KLOE Kl //3 [52] 


22.3 ± 10.5 


4.8 ± 5.2 


9.1 ±6.5 


-0.97 


+0.81 


-0.91 


Et 




KLOE Kl e3-//3 [52] 


25.6 ± 1.8 


1.5 ±0.8 


15.4 ±2.2 


-0.95 


+0.29 


-0.38 


average 


1 


KTeV Kl e3 [48] 


21.67 ± 1.99 


2.87 ±0.78 




-0.97 






tl 




KTeV Kl /x3 [48] 


17.03 ± 3.65 


4.43 ± 1.49 


12.81 ± 1.83 


-0.96 


+0.65 


-0.75 






KTeV Kl e3-;u3 [48] 


20.64 ± 1.75 


3.20 ±0.69 


13.72 ± 1.31 


-0.97 


+0.34 


-0.44 


average 


1 


NA48 Kl e3 [51] 


28.0 ±2.4 


0.4 ±0.9 




-0.88* 






^Eyi ilowl ^high) 


2 


NA48 Kl /u3 [54] 


20.5 ±3.3 


2.6 ± 1.3 


9.5 ± 1.4 


-0.96 


+0.63 


-0.73 






ISTRA+ K- e3 [50] 


24.85 ± 1.66 


1.92 ±0.94 




-0.95* 






(y,^)2Cfit 


3 


ISTRA+ K- fi3 [53] 


22.99 ± 6.42* 


2.29 ±2.29* 


17.11 ±2.25* 


-0.82* 


-0.12* 


-0.41* 


(y, ^)2Cfit 


4 



, I— I ^-^abliQC- kleasurements of Ki^ form-factor slopes. Values marked with an asterisk involve additional assumptions; see notes in text. 

□ □ □ □ 



